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1. Introduction 



The Large Hadron Collider (LHC) at CERN is now beginning operation, and will be the 
major source of data from the energy frontier for many years to come. The main goal of 
the LHC experimental program is to discover physics beyond the Standard Model (SM), 
as well as the mechanism of electroweak symmetry breaking. In the SM, this breaking is 
due to a single scalar field, the Higgs field |Q, ^, Similar fields exist in most extensions 
of the SM, such as the Minimal Supersymmetric Standard Model (MSSM). 

The discovery and the measurement of the Higgs sector will be a central piece of the 
experimental effort at the LHC. The production of the Higgs boson is dominated by gluon 
fusion through a top quark loop, for the whole relevant Higgs mass range [^. The next- 
to-leading-order (NLO) corrections to the gluon-fusion production cross section are very 
large, of the order of 100% ^, 0, The second most important contribution to the 
Higgs production cross section comes from the vector boson fusion (VBF) process, which 
proceeds at tree-level and receives much smaller QCD corrections [^, |10|, . 

Both of these channels participate in the phenomenologically interesting signal of pp 
H + 2 jets. The jets coming from the two processes have different angular distributions: 
well-separated and forward jets in the VBF case, in contrast to less separated jets and 
further central jet activity in the gluon fusion case. Thus, the cross-contamination can be 
reduced by imposing appropriate experimental cuts. A good theoretical understanding is 
also necessary to reduce the uncertainties coming from the backgrounds and interference 
effects, and to allow us to perform precision studies of the Higgs sector. 

For the gluon-fusion contribution, in order to make the computation more tractable, 
it is possible to make the approximation of a large top quark mass mt [12, 13, ^, |^. 



This approximation replaces the full one-loop coupling of the Higgs boson to gluons via 
a top quark loop, by an effective local operator C (mt) H G G^'^ ; thus it reduces the 
problem by one loop order. Formally, this approximation requires the Higgs mass to obey 
rriH <^ 2mt ~ 350 GeV. However, for inclusive Higgs production, the approximation works 
very well up to quite large Higgs masses, mn ~ 2mt or more, if G{mt) is taken to have 
the exact mt dependence from one loop |lj]. Here we are interested primarily in processes 
where the Higgs boson is relatively light, but because of the extra jet activity the partonic 
center-of-mass energy and final-state invariant masses may be large. It has been shown 
that the large-mt approximation is still valid for such configurations, as long as the jet 
transverse energies are smaller than mt jl^, [T6| . 

Several groups have computed various relevant quantities for the gluon-fusion contri- 



bution to pp ^ H + 2 jets, at LO, in the large mt limit |17, IT^ and with the exact mt 



dependence [20|; and at NLO accuracy in the large mt limit [21, The real corrections 
to this process, involving tree amplitudes for a Higgs plus five partons, were studied in [p3|[ . 
The interference between gluon fusion and VBF production has been computed, and found 



to be very small [24, 25[. Our results for the HqqQQ amplitudes can also be used to 



- 1 - 



calculate the one-loop interference between the color-singlet pieces of these two processes. 
We outline this calculation in section ^ 

New methods have been developed for computing one-loop amplitudes for multi-leg 
processes. Some are based on Feynman diagrams [26, |2^, ^ |3^, 33, 



35 



36] while others exploit generalized unitarity |37, 35, 41, 42, 43 1 and recursion 

relations [Q]. The methods based on Feynman diagrams have been employed to compute 
several quantities involving the Higgs boson to next-to-leading-order 



mm 



In the large rrit approximation, the NLO corrections to the production of a Higgs boson 
plus various numbers of jets at a hadron collider require one-loop amplitudes in QCD, with 
one insertion of the effective operator HGfj_yG^'^. We will refer to these amplitudes as 
one-loop Higgs plus multi-parton amplitudes. Because these virtual amplitudes contain 
infrared divergences, they are invariably presented using dimensional regularization (we 
take D = A — 2e), as a Laurent expansion in e through the finite 0{e^) terms. The 
complete set of such amplitudes for three external partons {ggg or qqg) was provided in 
ref. [E^]. Results for various numbers of legs and helicity configurations have appeared 



more recently [4S, 50, 51, 52]. In particular, the amplitudes with four gluons, all of the 



same helicity [49], and those with two negative and two positive helicities ]51, 52] have now 
been computed analytically, using techniques very similar to what we will employ here. 

The full analytic results for the one-loop Higgs plus four-parton amplitudes, for a 
complete set of parton helicities, have not yet appeared. However, the analytic expressions 
for the color- and helicity-summed NLO interference of one-loop and tree amplitudes have 
been presented for the four-quark case J^] , along with numerical results for the two-quark- 
two- gluon and four- gluon cases. These results have been incorporated into the NLO gluon 
fusion contributions to pp ^ H + 2 jets mentioned earlier ]^2|. In this paper we present 
analytic results at the amplitude level for the four-quark case, HqqQQ. We also give 
the two-quark- two- gluon amplitudes, with the restriction that the two gluons must have 
opposite helicity, Hqqg^g^ . 

In our method, the Higgs field H is rewritten as the sum of a complex field (p, and 
its complex conjugate i;^^. This has the advantage that the analytic structure of the two 
components is much simpler than in the total amplitude. Furthermore, we only need to 
compute the (/>-amplitudes because parity relates them to the (f)^ ones. Our technique 
for calculating these amplitudes consists of a unitarity-recursive bootstrap approach: by 
performing appropriate unitarity cuts we obtain all cut-containing terms of the amplitudes 
(logarithmic, polylogarithmic functions, and associated terms), while the rational terms 
are computed using on-shell recursion relations. In this process we only use on-shell lower- 
point amplitudes as input in our calculation. This greatly simplifies our task and allows 
us to obtain compact analytic answers efficiently. 

In slightly more detail, we employ quadruple cuts J^] to determine the coefficients 
of scalar box integrals in the amplitudes. The only scalar triangle integrals appearing 
in the amplitudes described here have one or two external massive legs, not three; the 
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coefficients of such integrals are fixed easily using the amplitudes' known infrared poles. 
The coefficients of bubble integrals are computed using the method of spinor integration 
via residue extraction 55]. As just mentioned, the rational terms are computed using 



on-shell recursion relations. Certain spurious poles arising in this technique are dealt with 



using an approach which is a hybrid between the cut-completion method |5^, 57, ^ and the 
method of evaluating the spurious pole residue of the cut part |59] which is implemented 
numerically in BlackHat Q. 

Our analytic expressions can easily be incorporated into one of the NLO computer 
programs for computing cross sections. They provide a fast evaluation of the amplitudes 
and are more stable compared to (semi) numerical approaches. In the computation of Higgs 
amplitudes with yet one more external parton (five in all), they can serve to check limiting 
cases, when two partons become coUinear, or one gluon becomes soft. In a numerical on- 
shell recursive approach such as BlackHat [^4|, they could provide a more important role: 
the four-parton amplitudes could be used as a fast analytic input for some of the terms in 
the on-shell recursion relations for the five-parton amplitudes. 

By performing the appropriate sum over colors and helicities of the interference between 
tree and one-loop amplitudes, we are able to confirm (numerically) the expression for the 
virtual part of the NLO cross section for HqqQQ presented in ref. (2^. At the same time, 
because our results are color decomposed, we can project them into a color-singlet channel 
with respect to the qq (and QQ) quantum numbers. The color-singlet channel can interfere 
with the electroweak VBF process [|60|. We have verified numerically that this interference, 



which is one of the two virtual contributions to the full interference, agrees with one [24| 



of the two recent computations of this (quite small) effect [ SI | . 



The color-singlet parts of the one-loop amplitudes may be of use in determining how 
frequently the gluon-fusion process produces events with large "rapidity gaps" that would 
survive typical central jet vetoes proposed to select the VBF process. Resummed estimates 
of the efficacy of such vetoes have been performed recently in ref. [^2| for example; however, 
there may also be important contributions at fixed order in a^- 

The paper is organized as follows: In section |2| we introduce some basic notions that 
simplify the task of computing the HqqQQ and Hqqgg amplitudes — the (p-cj)^ and color 
decomposition of Higgs amplitudes, and the spinor helicity formalism. In section ^ we 
outline the basics of the unitarity-recursive method at the core of the calculation. We also 
apply the technique to specific examples. In section |^ we present the full analytic answers, 
numerical results, and applications of our expressions. We describe checks that we have 
used to verify their correctness. In section ^ we present our conclusions and mention 
possible future directions. 



2. Notation 



In this section we introduce the basic notation used in the rest of this paper, as well 
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as some prerequisite notions. In particular, we decompose the Higgs amplitudes into (j) 
and (f)^ components with simpler analytic properties, describe the color decomposition of 
scattering amplitudes in terms of partial and primitive amplitudes, and recall the spinor- 
helicity formalism. 



2.1 The 



decomposition of Higgs amplitudes 



In the amplitudes we study in this paper, all external quarks are taken to be massless, and 
the Higgs couples to them through gluons. The coupling of the SM Higgs boson to gluons 
is dominated by an intermediate top quark loop, because the top is much heavier than the 
other quarks. In the limit of very large top mass, mt — > oo, the top quark can be integrated 



out, giving rise to the following effective interaction |12, 13 1 

C 



£int 
H 



(2.1) 



with the coefficient C given by C = as/{QT^v) = / {2A7t'^v), to leading order in q^. Here v 
is the vacuum expectation value of the Higgs field, v = 246 GeV. (The value of C is known 
to 0(af) |6^-) Our convention for the normalization of generators is tr T^T'' = 5°"^ ^ and 



Gpu = EaT^Gf,,, so that trG^, 



Tree-level amplitudes (not counting the top quark loop) for a Higgs boson plus multiple 



partons were ffi'st computed analytically for up to four partons in refs. [17, 18, 1£], and up 

by 



to five partons in ref. |64]. The structure of these amplitudes could be simplified ||6^ 
splitting the effective interaction Lagrangian into two parts, a holomorphic (self-dual) and 
an anti-holomorphic (anti-self-dual) part. In fact, certain "maximally helicity violating 
(MHV) rules" that had been developed for QCD tree amplitudes [|6^ could be applied 
straightforwardly to Higgs amplitudes after making this split |^, 67 1. Specifically, we 



consider the Higgs boson to be the real part of a complex field (p, with 
The interaction Lagrangian can then be rewritten as 



C 

2" 
C 



fiti G sD pv Gg]j + (j)^ trGAsofiu G^gjj 



(2.2) 

(2.3) 
(2.4) 



where the gluon field strength has been divided into a self-dual (SD) and an anti-self-dual 
(ASD) component, given by 



-(G^^ + *G'^^) 



'-^ASD 



*Qt,u ^ l^t^'^p-Gp^ . (2.5) 



From eq. ( |2.2[) and its conjugate we can reconstruct the scalar H and pseudoscalar A 
fields according to 



H = 6 + 



A 



1 



(2.6) 
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It follows from eq. ( |2.6D that the scattering amplitude for a scalar Higgs boson plus any 
number of partons can be obtained, at any loop order /, by the sum of the amplitudes with 
(f) and (j)"^, 

A^\h,...)=A'^\cP,...)+a1\<P\...), (2.7) 
with ". . ." indicating any arbitrary configuration of partons. 

Similarly, the amplitudes for a pseudoscalar A plus partons are given by 

1 



A^\a,...) = ^ A^\<P,...)-A^\cP\...) 



(2. 



recognizing that the constant C is different for the A case |50|. The relation between the 
(j) and (j)^ amplitudes is through parity, or complex conjugation of spinors. 



aHH^^, i^s 2^^ . . . ,n^") = (-1)"^-'' U«(<A, 1 



-h2 



,n 



(2. 



where the spinor products (ij) and [j i] are defined in eqs. ( 2.26| ) and ( 2.27| ), and 



n 



denotes the number of external antiquark-quark pairs |64]. In other words, to generate an 
amplitude with (j)^ from an amplitude with (j), one reverses the helicities of all quarks and 
gluons, and replaces all spinors (ij) with [ji]. Thus, it is possible to reconstruct the H 
and A helicity amplitudes by computing only their (/)-components, using parity to get the 
''^-components, and then assembling the two ingredients together. 

It is important to note that in the separation of any Higgs amplitude into a (p and 
a 0^ amplitude, all color and kinematic information {e.g., the momentum of the scalar 
particle) remains the same in the original and the component amplitudes. What separates 
is the self-duality properties of the amplitudes, and consequently their analytic structure, 
resulting in a simplification of the calculation. 



2.2 Color decomposition and color sums 



An important tool for QCD calculations is the color decomposition of amplitudes |68, 69 



7C]. It allows us to write down a color-ordered expression for any amplitude, which is a 
sum of products of color structures and uncolored functions of the kinematic variables, 
called partial amplitudes. The color and kinematic information is neatly separated in 
this way, and one has to compute only the partial amplitudes, which have simpler analytic 
properties. The partial amplitudes can be expressed in terms of yet simpler building blocks, 
called primitive amplitudes. Primitive amplitudes are color ordered; that is, they receive 
contributions only from planar one-loop Feynman diagrams with a fixed cyclic ordering 
of the external legs. Furthermore, a given primitive amplitude either contains a closed 
fermion loop (f), or it does not. If it does not, then the primitive amplitude is further 
characterized by how the external fermion lines are routed as they enter the loop: whether 
they turn left (L) or right (R) in the case of one fermion line; and according to leading- 
color (Ic) and subleading-color (sic) designations in the case of two separate fermion lines. 
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(c) 



Figure 1: Sample diagrams corresponding to the (a) leading color or "Ic", (b) subleading-color or 
"sic" and (c) fermion loop or "f" primitive amplitudes. The Higgs field can attach to any gluon 
line in these diagrams. 

Computing the primitive amplitudes will be the focus of this paper. First, we describe the 
decomposition of the (pqqQQ and (pqqgg amplitudes in terms of partial, and then primitive, 
amplitudes. 



2.2.1 (t>qqQQ 

Because the Higgs boson, and (j) and ip'^ fields, are color-neutral, they play no role in 
the color structure of the amplitude. The color decomposition of the HqqQQ or cpqqQQ 
amplitude is identical to the decomposition of the four-quark qqQQ amplitude ||7l|. There 
are two independent color structures, corresponding to the two independent ways we can 
contract the color indices of the quarks. At tree level, the coefficients of the two color 
structures are simply related. 



Af\(^, Ig, 2„ 3q, 4q) = Cg' A'^>{(t>, 1,-, 2„ 3q, 4, 



,2 /.(O), 



where the tree helicity amplitudes A^^\(p,lq,'^qT'iQ,4:Q) are given in eq. (A^), and Nc 
denotes the number of colors, = 3 in QCD. 



(2.10) 



At one loop, the color decomposition is 



^^cCC^4;l(<^, 15,2„3q,4q) 



+ 5J^C^4;2(</',1,-,2,,3q,4, 



where we have also extracted a factor from the loop integrals, 

1 r(l + e)r2(l-e) 



cr 



(47r 



,2-e 



r(l - 2e) 



(2.11) 



(2.12) 



Formulas ( 2.10 ) and ( |2.1lD apply to the case of different quark flavors, q Q- The am- 
plitude for identical flavors, q = Q, is found from the unequal- flavor formula by subtracting 
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the same formula with the labels for q and Q exchanged, 



A. 



(0/ 



l5,2g,35,4g) =4')(</),1,-,2„3q,4 



A, 



(0/ 



Igi 4g, 3q, 2c 



(2.13) 



Note that the helicities of q and Q must be the same in order to get a nonvanishing 
exchange term. Of course if there are two identical quarks in the final state, there is also 
an identical-particle factor of ^ in the phase-space measure. 

The partial amplitudes ^4;! and ^44. 2 can in turn be expressed in terms of primitive 
amplitudes, using the results of refs. |71| for the analogous amplitudes, e"'"e~ — > qqQQ. 
Because the e~^e~ pair and the scalar (p are both colorless, the color structure is identical 
to our case. The results could be given in a helicity-independent form; however, we list 
the two different helicity cases separately, so that we can use relations among the primitive 
amplitudes in order to minimize the number that appear: 



A4;i(</., l5,2+,3t 4 



1 r 2+ .3+ 4" 



(2.14) 



iV2 



2 r 



4^(</),ir,2+ 3+ 43)+4^(</>, ir,2+ 4;i,3^ 



4 

sic 



n 



,lr,2+,3+,45) + ^4(^,lr 2+ 3+4 



^4;i(</',l,-,2+,3^,4+) = 4^(</<,lr,2+,3^,4+) 



(2.15) 



4^(</., ir , 2+ 3;i, 4+) + 4^(0, 2+ 4+ 3, 



+ -^"^ti^, 1," , 2+, 4+, 3q) + ^4(<A, 1,- , 2+, 3^, 4+) 



^2(0, l,-,2+,3t, 45) 



A, 
+ 



4 (ri -^q ) -5 

1 



Af2 



1 r ?+ 4T .s 

4^(</.,ir,2+, 4^,3+) +4^(0, 1^,2+, 3+ 4^ 



(2.16) 



+ ^AticP, ir , 2+ 3t , 45) - ^LA.i^, ir , 2+ 3t 45) 



/ /if, 



^4;2(0,l5,2+,3^,4+) = Alf(0,l^r,2+,4t,3Q) 



(2.17) 



+ 



4^(,/., ir, 2+, 4+ 3q) + 4^(0, ir, 2+, 33, 4+) 



Ar2 

c 

]V2^4''('/'> Ig . 2+, 4t , 3q) - jf 



4(<A,l„-,2+ 3;:,4+). 



Here , Af'^, and ^4 describe respectively the leading-color, subleading-color and fermion- 
loop primitive amplitudes. The number of massless quark flavors is denoted by nj. The 
quantity ^4 (<;/>, 1^, 2+, 4^, 3q) refers to the primitive amplitude A4 (</>, 1^, 2+, 3^, 4^) with 
the labels on legs 3 and 4 exchanged. Sample Feynman diagrams corresponding to these 
amplitudes are shown in figure ||. 
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1 4 1 4 1 4 

(a) (6) (c) 




Figure 2: Sample diagrams corresponding to the (a) qqgg L, (b) qqgg R, (c) qgqg L, (d) qqgg 
fermion loop, and (e) qgqg fermion loop primitive amplitudes. The Higgs field can attach to any 
gluon line in these diagrams. 



The virtual part of the unpolarized NLO cross section for HqqQQ requires the sum 
over both helicities and colors of the interference of the tree and one-loop amplitudes. The 
color sum is straightforward to work out from eqs. and ( ^ITD , 



colors 



= 2C^crgHN^-^)NcRe[Af'>*{H,lg,2„3QAQ)A4-i{H,^,2„3QAQ) ■ 

(2.18) 



The same formula also applies with H replaced by A. The formula for the case of identical 
quarks q = Q follows from eq. ( ^.13| ). 



2.2.2 cl)qqgg 

Similarly, the color decomposition for the (pqqgg amplitudes is identical to that for the 
process e~^e~ — > qqgg |72|. The tree amplitude is given by 



Af\^,lg,2„3A) = Cg^ E {T''<'^T''^''^)^l' Af\^A^,2„a{3),a{A)) , (2.19) 



where the tree helicity amplitudes ^4'^'* (0, Ig, 2g, 3, 4) are given in eq. ( [A.5| ) 
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The one-loop amplitude is decomposed as 



4'^ (0,15,2,, 3, 4) = Cg^cr 



(TM3)r«'^(4)).nyl4.^(^,i_,2„a(3),a(4)) 
+ <5«3a4A n^^.3(^^l_^2,;3,4) . (2.20) 



The partial amplitudes ^4-1 and ^453 are given, in a helicity-independent fashion, by 

1 



^4;i(</', I5, 2g, 3, 4) = (</., 1,-, 2g, 3, 4) - ^vlf (</>, 1,-, 2,, 3, 4) 



+ ^4(,/.,l,-,2„3,4) 



(2.21) 



A4;3(</', I5, 2g; 3, 4) = Al{(t>, I5, 2„ 3, 4) + (<^, I5, 2,, 4, 3) + Al{ct>, I5, 3, 2„ 4) 

+ ((/>, 1,-, 4, 2„ 3) + ((A, 1,-, 2„ 3, 4) + .4f (0, 15, 2,, 4, 3) . 

(2.22) 

Here we choose to label the leading- and subleading-color primitive amplitudes by "L" 
and "i?" (corresponding to fermion lines turning left or right upon entering the loop) 
to be compatible with the notation used in ref. |50|. Again ^4 denotes a fermion- loop 
contribution. Sample Feynman diagrams corresponding to these primitive amplitudes are 
shown in figure |^. 

The virtual part of the unpolarized NLO cross section for Hqqgg requires the sum over 
both helicities and colors of the interference of the tree and one-loop amplitudes. The color 
sum can be expressed in the same form as that for e+e~ — > qqgg f?^. 



NLO 



colors 



2C7^cr/(A^^ - 1) Re<j *(ii', I5, 2„ 3, 4) {N^ - l)^4;i(^, I5, 2,, 3, 4) 



(2.23) 



- Ai,i{H, Ig, 2g, 4, 3) + A^-siH, I5, 2,; 3, 4) 

+ {3^4}. 

The same formula also applies with H replaced by A. 



2.3 Spinor helicity formalism 

Primitive amplitudes depend only on kinematic variables. They are functions of the exter- 
nal momenta of the Higgs boson, = kn, and of the four partons, fej, i = 1, . . . , 4. These 
momenta are all outgoing, by convention, so momentum conservation and the on-shell 
conditions read. 



k^ + y^kj = , 

i=l 

T 2 _ T 2 _ 2 



kf = 0. 



(2.24) 
(2.25) 
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A very convenient representation of the amplitudes is in terms of spinor inner products, as 
reviewed e.g. in refs. |6^, |7^. Let u±{ki) be a massless Weyl spinor of momentum ki and 
positive or negative chirality. The corresponding two-component spinors are often denoted 
A" and Xf. The spinor products are defined by 

{ij) = XfXja = (r = u-{ki}u+{kj) , (2.26) 
[ij] = ~Xia~Xf = {i^lD = u+{k.i)u-{kj) . (2.27) 

We use the convention [ij] = sgn(A:°A:°) {j i)* , so that 

{ij)[ji] = 2ki-kj = Sij. (2.28) 

For real momenta, and up to a complex phase, the spinor inner products are square roots 
of the corresponding kinematic invariant Sij = {ki + kj)'^. 

Three-parton invariant masses are defined by, 

Siji = ih + kj+kif = {i j) [j i] + {j I) [I j] + {i I) [I i] . (2.29) 

We also define the spinor strings, 

{a\i\b] = {ai)[ib] , {a\ii + j)\b] = {ai)[ib] + {aj)[jb]. (2.30) 

Strings involving the Higgs momentum, such as {a~\kH\b~) , could also show up; however, 
they can be eliminated in favor of strings such as in eq. ( 2.30| ) using momentum conserva- 
tion. 

Besides momentum conservation, the other two spinor product identities used to sim- 
plify expressions are antisymmetry, 

{3i) = -{ij), [ji] = -[ij], (2.31) 

and the Schouten identity, 

(a 6) {cd) = {ad) {cb) + {ac) {bd) , (2.32) 

[ab] [cd] = [ad\ [cb] + [ac] [bd] . (2.33) 

Composite spinors can appear in these products as well and they can be handled similarly, 

{a]P^.Jb] {cd) = {ad) {c]P^.Jb] - {ac) {d]Pi.Jb] , (2.34) 
{a]P,,„j]b] [cd] = {a]Pi,„j]d\ [cb] + {a[Pi,„j[c] [bd\ , (2.35) 
{a[P^...j[b]{c[Pi.Jd] = {a[Pi...j[d]{c[Pi...j[b] - Ptj {ac) [bd] , (2.36) 

where Pi...j denotes an arbitrary momentum sum, Pf^ ■ = X^^^j km- The spinor products 
defined above form the building blocks for the primitive amplitudes and the basis for our 
calculations in the next sections. 
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3. Unitarity and recursive bootstrap method 



The method we employ in our paper combines the (generalized) unitarity method |73, 74 



72, 53, 54, 55] with on-shell recursion relations [^] operating at one loop |76, 77 



Using general integral reduction formulae ||7^, any dimensionally-regulated one-loop 
amplitude An^ with massless internal lines can be decomposed as, 

^« = Cn + Rn, (3.1) 

where the cut part Cn is a linear combination of scalar integrals |7^ , 

Cn = Y.diii+Y,ci4+J2^iA, (3.2) 

i i i 

and Rn denotes the rational part, which contains no branch cuts (in four dimensions). 
The box integrals II, triangle integrals Xg and bubble integrals T2 are well known (see 
for example refs. Q). Hence determining C„ is equivalent to computing the 

coefficients of the respective integrals, di, Ci and hi. These coefficients are found by taking 
(generalized) four-dimensional unitarity cuts of the amplitude in various channels. The 
rational part Rn will be calculated by utilizing on-shell recursion relations, which requires 
only information about lower-point loop amplitudes and tree amplitudes. (For a review of 
the on-shell approach, see ref. 



In the remainder of this section, we outline the various ingredients in the method, and 
apply them to specific primitive amplitudes for (pqqQQ and 4>qqgg. The results for these 
amplitudes are collected in section |^. 

3.1 Generalized unitarity for box coefficients 

In this paper, the generalized unitarity method, along with complex spinor integration |53 



|5^ , pq ] is employed in a twofold manner for the calculation of box functions and bubble 
(single log) functions respectively. 



In ref. |53] it was shown that using quadruple cuts, it is possible to reduce the task of 
calculating the coefficient of any box function into that of calculating a product of four tree 
amplitudes. This method can be applied to box functions with any number of external 
masses by taking advantage of the generally non-vanishing behavior of the all-massless 
three-point amplitudes when momenta are complexified. A three-point amplitude with 
one leg carrying opposite helicity from the other two can have either an MHV or an anti- 
MHV (sometimes also denoted MHV) representation while satisfying energy-momentum 
conservation. The complex momenta allow for all four loop momenta of a box diagram to 
become on shell, freezing the cut integral completely and simplifying it into an algebraic 
product of tree amplitudes. 

For example, the coefficient of the box integral in figure ^ is given by 

d = —V f dHi 5^+\tl)5^+\ll)6^+\il)5^+\ll) Af^A^^^AfAf , (3.3) 
Als2-4 J 
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Figure 3: Evaluation of a box function coefficient using quadruple cuts. The four on-shell condi- 
tions on the loop momentum reduce the integration to an algebraic product of four tree amplitudes. 

where £2 = — Ki, £3 = ii — Ki — K2, £4 = £1 + K4, and each Af'' corresponds to the tree 
amplitude in the respective corner with massive (or massless) momentum Ki (i = 1 . . . 4) . 
The quadruple cut (or leading singularity) of the scalar box integral is given by 

AlsX4 = I d% sM(£l)6M(£l)6M(il)sMi£l), (3.4) 

from which we obtain 

where the summation is over the two discrete solutions ^J, a = 1,2, of the loop- momentum 
localization constraints, 

£l = £l = £l = £l = 0, (3.6) 
and over all possible helicities h of internal particles propagating in the loop. 



3.1.1 A two-mass box coefficient 

Let's consider, for example, the leading-color primitive amplitude ^4 ((/>, 1^, 2+, 3^, 4q). 
One of the box coefficients that needs to be determined for this primitive amplitude is that 
of the "easy two mass" box integral, with diagonally opposite massive legs having mass 
and S23- Because all the two-mass boxes have m'jj as one of the two masses, we denote 
this coefficient, using the other mass, as d'^^'^- This box integral is defined by the clustering 
of the five external particles into the four legs of the box: (1) (23) (4) ((/>). The associated 
quadruple cut is shown in figure |^. 

The tree amplitude with <p and two gluons vanishes unless both helicities are nega- 
tive [S5|. It is easy to see that this fact, plus fermion helicity conservation, forces the 
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unique assignment of intermediate helicities shown in figure |^. The three-point vertices 
containing legs 1 and 4 are only nonvanishing for one of the two solutions to the quadruple 
cut conditions (p^). Equation (|3.5D then becomes, 



(3.7) 

Using the expressions for the tree amplitudes, we get 

rJ^^^-'\(/f(f')\'\><\tllll£^( [23]^ \ [i'A-U)]' 



The positive chirality spinors for the three-point vertices containing legs 1 and 4 are 
proportional, 

A^^ oc A^2 oc Ai , Xi^ oc Xg'^ oc A4 . (3.9) 

Using this fact, along with momentum conservation relations, we can eliminate all explicit 
loop momenta from eq. (^J 



23 2^ ^ [1^2] [^2 2] [34] [^4 4] 
^ 1 2 ([^2l](14)[44])2 



2' ' [14] [42] [34] [44] 

1,. ,,2r.^l2 [1^2](44) (14) [44] 

2^ ^ ^ ^ [2 4] (44) (14) [4 3] 
l,.,^2rg ..12 [11(2 + 3)14) (1|(2 + 3)|4] 
2^ ' ^ ^ [23] (34) (12) [23] 

^(4|(2 + 3)|1](1|(2 + 3)|4]^^^^. (3.10) 
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Figure 5: Quadruple cut for the evaluation of the one-mass box coefficient d\™^ appearing in 
Alf(0,f,r,2+,3t,4Q). 

In section ^ we express the result, not in terms of the scalar box integral Xl™*^, but in 
terms of the infrared-finite box function Ls?_™^ (si23, S234; S23, "t-//) • These are related by 
eq. ( |B.6|) . After removing a factor of cp associated with eq. (|2.11| ), and using the identity 

S123S234 - S23m% = (4| (2 + 3) 1 1] (1 1 (2 + 3) |4] , (3.11) 

the coefficient of the box function Ls?.™^ (•S123, 5234; S23) "i^) in A^^{(l), 1^, 2+, 3^, 4q) is 

Dir ^ — = i Jif.,, = -4°^(<^,i,^,2+,3t,4Q). (3.12) 

S123S23A - S23mjj (12) (3 4) v y y >^ 

It is a general feature of every primitive amplitude presented in section |^ that all easy- 
two-mass box functions have coefficients equal to the negative of the corresponding tree 
amplitude. We therefore collect all the Ls?.™*^ functions into 'V functions, which also 
contain the infrared and ultraviolet poles in e, since the latter have to be proportional to 
the tree amplitude as well. 

3.1.2 Calculation of a one-mass box coefficient 

As our second box example, we compute for the same primitive amplitude the coeffi- 
cient of the one-mass box function Ls_i (si2, S23; S123)! associated with the external leg 
clustering (1)(2)(3)(4(/)). The quadruple cut is depicted in figure We label the box 
integral coefficient by ^123- Again the vanishing of a tree amplitude involving cj), in this 
case A'^\(/),£f , —£~^q,4q), along with fermion helicity conservation and the vanishing of 

A-^4\q~ i 9^ J 9^)^ forces the unique helicity assignment shown in the figure. The cut 
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4 - 

Figure 6: Sample (j>qqQQ triple cut, illustrating the vanishing coefficient for a three- mass triangle 
integral, caused by the tree amplitude in the lower right corner associated with the external invariant 
S34. Consequently, hard- two-mass box functions are absent from the amplitude as well. 



evaluates to 

"2^1^ (Fm) ^^ [ {{-i2) 2) ) ■ ^ ^ ^ 

After some spinor product manipulations similar to eq. ( p.lO ), we have 

which yields for the coefficient of the function Ls_i (si2, S23; ■5123), 

512^23^''' '(12) (34) vr,^5,-,,"Q 



Du3 = -^<3 = ^,lwLv = -4°^(0,l5,2+ 3+ 4.). (3.15) 



The result is again proportional to the tree, up to a sign. This property holds for all 
the one-mass box coefficients Ls_i in the (pqqQQ primitive amplitudes, apart from the 

leading-color piece of the ( — I h) helicity configuration. However, in the (pqqgg primitive 

amplitudes it is typically not true. 



3.2 Absence of three-mass triangles 



The one-loop (f) amplitudes that we consider in this paper have the property that the triple 
cuts associated with three-mass triangle integrals all vanish. This general feature holds 
because the tree-level four-parton amplitudes in QCD, A^'^^ (1, 2, 3, 4), and those with an 
additional 6 boson, Af^((/>, 1,2,3,4), each require two negative helicities to be nonvanishing. 



A nonvanishing product of three such amphtudes, as required for a triple cut, imphes six 
negative hehcities. Three of the negative hehcities are associated with the three cut hues, 
so there must be three external negative hehcities. However, the one-loop (j) amplitudes 
we consider here only have two negative hehcities. An example of a vanishing triple cut 
is shown in figure ^. The same argument implies that all "hard-two-mass" boxes, with 
two adjacent massive legs, must vanish: For any hard- two-mass box quadruple cut, one 
can remove the cut between the two adjacent massless legs, relaxing the quadruple cut 
into a (vanishing) triple cut. Because the (/> -|- 4 parton amplitudes obviously contain no 
three-mass or four-mass box integrals, only one-mass and easy-two-mass box coefficients 
have to be computed here. 



The amplitudes (j)g~ g~ |^T| and (pg^ g^ g^ g^ [£2[ also contain only one-mass and 
easy-two-mass boxes, and no three-mass triangles, for the same reason, insufficiently many 
negative hehcities in the triple cuts. Interestingly, the amplitude (jyg g g g^ computed in 



ref. |49| also contains only one-mass and easy-two-mass boxes, and no three-mass triangles, 
for the opposite reason, a paucity of positive helicity gluons. On the other hand, the 
primitive amplitudes for 4'q~q^g~g~ and (j)g^ g~ g~ g~ , which have not yet been computed 
analytically, will contain three-mass triangles and hard-two-mass boxes. 

For the coefficients of two-mass and one-mass triangles, the above triple-cut vanishing 
argument does not hold. The existence of a massless external leg implies that one of the 
tree amplitudes can be an MHV three-point amplitude, which contains only one negative 
helicity, not two. However, the two-mass and one-mass triangle integrals, given in eqs. (|B.7| ) 
and ( |B.8| ), contain single log terms at order 1/e. Because of this, their coefficients are 
completely determined by the known infrared poles of the amplitude, so they do not have 
to be computed separately. The remainder of the work to compute the cut part of the 
amplitude involves determining the coefficients of bubble integrals. 



3.3 Unitarity and spinor integration for bubbles 

In the case of the ordinary two-particle cuts used to determine bubble coefficients, there 
are not enough constraints to fully localize the cut integral. The cut contains a residual 
phase-space integral. The method we use in this case was proposed in refs. [54, 55 1 and 
consists of writing the cut loop momentum integral as an integral over spinor variables 
£ = \i) and £ = \£] [^^. The integrand can be transformed into a total derivative in £, 
leaving a single integral over £ which can be evaluated by residue extraction. 

For the cut shown in figure ^, the coefficient of the logarithm will be given by an 
integral of the form 

/ dHi6^+\£l) 6^+\el) Af^Af^ = dtt {i di) [£d£] 6 {P^-t{£\P\£])f {£,£), (3.16) 



where £ can be either £i or £2, whichever is more convenient. The function f{£, £) represents 




Figure 7: Example of the evaluation of a bubble (single log) function using an ordinary two-particle 
cut. Integration over the loop momentum is required in this case. 



the product of the two tree amplitudes and is in general a sum of terms of the form 

, with Qfc/P, (3.17) 



After performing the integration over t and partial fractioning using Schouten identities, 
we can always bring the remaining integrand into a form where we can take advantage of 
the identity 

and convert it into a total derivative with respect to (In the special case of n = 0, 
the spinor \r]] appears only on the right-hand side of eq. ( |3.18| ); hence it can be chosen 
arbitrarily.) Then we can evaluate the integral over \£) by calculating the residue for each 
pole. The case of multiple poles was also examined in ref. [p5|. 



In this process one also detects the coefficients of box integrals sharing the same cut. 
These are the terms that scale like 1/(^|P|£] after partial fractioning of the integrand. They 
can serve as an independent check of the box coefficients determined by the quadruple cuts. 



3.3.1 Calculation of a bubble coefficient 



As an example, we compute the coefficient of the bubble integral X2(si23), or equivalently, 
of the single logarithm ln(— S123), in the primitive amplitude Ajf (</), 1^, 2+, 3^, 4g). The 




Figure 8: Two-particle cut for the evaluation of the bubble function coefficient 6123 of 

^5' 



A!f((/),ir,2+,3+ 4, 



corresponding two-particle cut is shown in figure and the cut integral is given by 



i3/dLIPS -M: X (1(-«>^P3) 



/ 



(£2 4) ; (1 2) (2 3) (3 i-h)) ii-h) i-ii)) ii-ii) 1) 
dLIPS 1)^(^1 4)^ 



Here dLIPS stands for the Lorentz-invariant phase space measure. 



By multiplying both the numerator and denominator of the integrand in eq. ( 3.19| ) 



with factors of [£2 £1] we can eliminate one of the loop momenta (in this case we choose to 
eliminate £1), leaving us with an expression that depends only on the other loop momentum 
(in this case, £2)- Furthermore, we rewrite the integral over the Lorentz-invariant phase 
space as an integral over spinor variables, as in eq. ( p.l6| ), 

dLIPS = j tdtj {£2 M2) [£2 M2] 5{P^-t {£2\P\£2]) , (3.20) 

with P = P123 = ki + k2 + ki,. We also need to track factors of ^/t from rescaling the £2 
spinors in the integrand. Performing these steps yields 

fit ! i.l.M,)V,M,\8(P^-tmP\l,\) (3.21) 







We can readily integrate over t, eliminating the (5-function, to get 



/ 



i.kim.M,\ p^i'-^'m^ . (3.22) 
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The two-particle cut we have considered detects not only bubbles, but also boxes that 
have cuts in this channel. By using Schouten identities we can rearrange terms so that we 
separate the bubble contributions from the box contributions, 



/ 



(£2 M2) \h dh] 



P2(£2l)(l|p|^2](^2 4) 2P2(^2 1)(14) 



+ \ (3 23) 

^ (^2|P|^2](12)(^2 3)(1|P|^2](^2 4) f ■ ^ ■ ^ 



We identify the box contribution as the last term in eq. (|3^ ), with the 1/ (^2|P|^2] depen- 
dence. Because the box contributions are obtained straightforwardly from the quadruple 
cuts, we can safely discard them (or use them as an independent check, but we won't do 
so here). The remaining part will be the coefficient of the bubble 2r2(si23), given by 

Transforming the integral over \t2\ into a total derivative using the general formula (|3.1 j ), 
with n = 0, 1 for the two terms in eq. ( 3.24 ), we obtain 



'123 



(£2 dl2) W2 a,,] - ^^^1^1^^^, + ^^^i^i^^j 3^ ^^^,^,3^ . (3.25) 



In this last step we have chosen the value \rf\ = |3] for the arbitrary spinor I?;] appearing in 
the n = term after its transformation into a total derivative. 

At this point we are ready to evaluate the integral over 1^2) by computing the residues 
of the poles of the integrand. We only have simple poles occurring for \l2) = |3) in the 
first term, and for 1^2) = -P|3] in the second term of eq. ( 3.25| ). Note that our choice for 
\rf\ eliminates a pole for \l2) = |3) in the second term. After substituting and simplify- 
ing, using (3|P|3] = S123 — S12, we get for the coefficient of the single log ln(— S123) in 



-B123 = -ib 



123 



1(1 2) [23]^ (34) ^[23]{IA) 



2 (Si23 - Sl2)^ S123 - S12 



(3.26) 



which is a rather simple final answer. 



The procedure outlined above for the computation of P123 is a typical example of the 
steps that have to be carried out to calculate any bubble function coefficient. It has been 
automated and implemented in Maple and yields a fast analytical evaluation of these cuts. 



3.4 On-shell recursion 



The (four-dimensional) unitarity technique can give us the cut-containing parts of the 
amplitudes (terms associated with functions with an imaginary part), but not the parts that 



are rational functions of the kinematic variables. However, these terms can be determined 
from their analytic properties as well, namely their factorization poles. On-shell recursion 



relations, developed first at tree level |84, ^] and later at one loop 77, 57, 55, |85| , 
exploit the known factorization behavior and have greatly simplified the task of calculating 
these terms. 



At tree level, all amplitudes are rational functions and one can consider |7^ a complex 
shift of any two of the external momenta j and / of an amplitude An , given by 

Xj Xj - zXi , Xi ^ Xi + zXj . (3.27) 

This [j, I) shift preserves momentum conservation as well as the massless conditions for the 
momenta kj and ki, which are now modified as follows, 

^ - iohno , - K + |(ii7no • (3.28) 

The shifted amplitude An{z) is an analytic function of z with only simple poles, which are 
associated with factorizations of An{z) onto lower-point tree amplitudes. 

Provided that An{z) — > as z — > oo, the integral of An{z)/z over the contour at infinity 
vanishes. This integral is also given by the sum over the residues at the poles for finite z. 
Therefore, the unshifted physical amplitude we wish to compute is given by the residue at 
z = 0, 

An = An{0) = - V Res^^. (3.29) 

Z = Za Z 

poles a 



Each pole Za in eq. ( 3.2S| ) is associated with a physical factorization channel of the ampli- 



tude. Factorization allows the evaluation of the residue, leading to the tree-level on-shell 
recursion relation |84, 75 1 

= E E Ai{Zrs)^A-^\z.rs) , (3.30) 

where h = ±1 denotes the helicity of the intermediate state carrying momentum Pr...s- The 
double sum over r and s is over partitions of the external legs into two sets (contiguous 
with respect to the color-ordering), for which the shifted legs j and / lie on opposite sides 
of the pole (j E L and I £ R). In the case of (p amplitudes, because (p is uncolored (and we 
do not shift the (p leg), it can appear on either the L or R side. The tree amplitudes on 



each side are evaluated at the complex momenta ( 3.28 ), shifted hy z = Zrs, where 



K s 

= {j-\Prl\l-) ^^-^^^ 

is the solution to the condition g(z.rs) = 0. 

At the one-loop level the situation is more intricate. Figure |^(a) shows schemati- 
cally the pole structure of a typical tree amplitude An\z). As shown in figure ^(b), the 
shifted one-loop amplitude An\z) can in general have not just simple physical poles, but 
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Figure 9: Analytic structure in the z plane of (a) tree amplitudes ^n\z), (b) one-loop amplitudes 
(z), and (c) the rational part of one- loop amplitudes Rn{z). The contour of integration is shown. 
The (green) dots in (a) and (b) represent physical poles. The additional (red) dots displayed in (c) 
represent spurious poles. 



also branch cuts in the complex plane, as well as double poles \7t]. In addition, it more 
frequently has non- vanishing behavior at infinity. 

Branch cuts would result in the need to evaluate discontinuities along an integration 
contour, as shown in figure |9|( 



_ To avoid this, one can use the decomposition (3J) of the 
one-loop amplitude An'' into a cut part C„ and a rational part and then work with the 
shifted rational part Rn{z), instead of analyzing the behavior of the full shifted one-loop 
amplitude An\z). Although Rn{z) is a rational function of z, only containing poles, some 
of these poles are spurious. The spurious poles are represented by the additional (red) dots 
in figure ^(c). Unlike the physical poles, the spurious poles are not associated with physical 
factorization channels. Their contribution to the entire shifted amplitude An\z) cancels 
between the shifted cut part Cn{z) and the shifted rational part Rn{z). 



^ in 



One type of spurious pole arises from the existence of terms such as (lnr)/(l — r) 
Cn- Here r is a ratio of two kinematic invariants that differ by a single massless leg. For 
example, eq. (|3.26| ) displays a factor of (si23 — su)^ in the denominator of the coefficient 
of ln(— S123) in the primitive amplitude ^4 ((/>, 1^, 2+, 3^, 4q). It corresponds to a term 
(lnr)/(l — r)^ with r = 5123/512- Under many choices of shift ( 3.27 ), r will become a 
nontrivial function of z. In this case a spurious pole, located at the solution to r{z) = 1, 
will be generated for Cn{z), and a compensating one for Rn{z). 

An analytic method for handling the contributions of spurious poles was developed in 
a number of papers |56, 57, 58, 85 1. The method has also been applied to Higgs boson am- 



plitudes 1 51, 52 1 . It consists of the following approach: We assume that the cut-containing 
pieces Cn have been obtained using methods such as those described in the previous sub- 
section. Then, for a general shifted one-loop amplitude we can write 



(Z) = Cn{z)+Rn{z). 



(3.32) 
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Our goal is to compute the rational terms We can absorb spurious singularities present 
in Rn into the cut-containing pieces by rewriting 

A^^\z) = Cn{z) + Rn{z) , (3.33) 

where the completed- cut terms Cn{z) are free of spurious singularities, as are the remaining 
rational terms Rn{z). 

To absorb all spurious singularities located at solutions to r(z) = 1, we make substi- 
tutions in Cn of the form 

In r In r -|- 1 — r , , , , 
2 - n_^,2 ^ Li(r), (3.34) 



(l-r)2 (1-r) 

Inr Inr — (r — l/r)/2 

(l-r)3 ^ (l-r)3 



= L2(r), (3.35) 



where r represents the ratio of any two kinematic invariants that differ by a single massless 
leg. The amount by which the completed-cut terms have changed in this process is given 
by the rational completed-cut terms, 

CRn{z) = Cniz) - Cniz) . (3.36) 

Consequently, 

Rn{z) = Rn{z)-CRn{z). (3.37) 

One can then consider the contour integral at infinity for Rn{z)/z. 

Provided that all spurious poles are removed from Rn{z) by the substitutions ( 3.34| ) 
and ( |3.35|) , this contour integral leads to an equation analogous to the tree-level recursion 
relation ( ^.29)) , featuring residues only at physical poles, 

^„ = ^„(0) = - y Res^^^^. (3.38) 

poles a 

These residues can be split into two sets of terms, using eq. ( p.37| ). The first set consists of 
the recursive diagrams, associated with residues of Rn{z); it can be evaluated analogously 
to the tree- level recursive diagrams ( |3.30 ): 

R^^ - y Res ?^ 

poles a 

— 1 R{kr, . . . , kj, . . . ,ks, —Pr...s) ^^'^^^s+lj ■ ■ ■ , kl, . . . , kr-l, Pr...s) 



h r,s 



+ vl*-'^-' {kr, . . . ,kj, . . . ,ks, Pr.^s) ■ ■ ■ ,kl, . . . , kr-l, Pr...s) 



h.iPAPl.s 

(3.39) 



A^^"^ (kr, . . . ,kj, ■ ■ ■ ,ks, Pr...s) p2 A^^\ks-\-l, . . . , kl, ■ ■ ■ , kr~l, P^ g) \ ■ 



r...s 
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The recursive diagrams are computed from the rational parts R of lower-point loop am- 
plitudes, and lower-point tree amplitudes A^^\ as well as the rational part of the factor- 
ization function which only enters for multi-particle poles [^6| (and not for collinear, 
two-particle channels). Just as at tree level, (f) can appear on either side of the pole. 

The second contribution from the physical poles consists of the overlap terms, 

On= y Res^l^^. (3.40) 

poles a 



They correct for the difference between Rn{z) and Rn{z) in eq. ( 3.37] ). In section ^.4.2 we 



will describe a modification of this procedure that can be used when the cut-completion 
described above fails to remove all spurious poles. 

Finally, we have to consider the potential contributions to the integral from infinity, 
because Rn{z) = An\z) — Cn{z) may not vanish as z ^ oo. In the case of the HqqQQ 
amplitudes there is always a shift that ensures a vanishing behavior of Rn{z) for large z, but 
there is no guarantee that this is always the case. In fact, some of the shifts used to compute 
the Hqqgg amplitudes have non-vanishing large z behavior. Usually it is straightforward to 
compute the z — > cxd limit of An{z) and Cn{z), denoted by Inf^„ and Inf C„ respectively. 
In some cases, a pair of shifts is necessary [^] (the original shift plus an auxiliary one), but 
that was not required for the amplitudes computed here. Putting together all the pieces, 
the full answer is given by 

Ai-i°°P = Cn + R^ + On- Inf Cn + Inf^n • (3.41) 



3.4.1 Calculation of rational parts for A^^{4i,l^ ,2+,3^,4q) 

To illustrate the calculation of the rational parts, we consider the primitive amplitude 
^4 ((/>, 1^, 2+, 3^, 4q). After completing the cut terms to form CR4, the first step is to 
choose a pair of legs [j,l) to shift according to eq. ( |3.27| ). Then we compute the recursive 
diagrams R^ and overlap terms O4, as well as any contributions from infinity (Inf terms) 
under this shift. 

For the case of ^4 (<;/>, 1^, 2+, 3^, 4^), the only parts of the cut terms that need com- 
pleting, to remove spurious singularities, are the single log terms. From the first term in 



eq. ( 3.26 ) we see that the function Li( Jg^^ ), as defined in eq. ( 3.34 ), should be introduced 



to remove the singularity as S123 — > si2- Similarly, the coefficient of ln(— S234) (which is 
related by symmetry to that of ln(— S123)) requires the function Li(^^^). These functions 
are collected in eq. ([4.^). From the rational parts of the Li functions we obtain CR4, 

CR - ' (3 4) [2 3]^ i (12) [2 3]^ 

^ 2 [1 2] (si2 - S123) 2 [3 4] (S34 - S234) ■ ^ ' ' 



-23- 



2 + 1 



2 + 1 





3 + 




3 + 



Figure 10: Diagrams needed to evaluate the recursive diagrams R^' of A]^{(j), 1^ , 2+, 3^, 4q). 

At this point there are no spurious poles left in C4, so we can proceed to choose a 
complex momentum shift ( |3.27D . For (1^, 2|j", 3^, 4q) we choose the [4,2) shift, namely 



A4 — > A4 



A2 ^ A2 + zX4^ , 



(3.43) 



or equivalently, 



|4] = |4]-z|2] 



z|4) . 



(3.44) 



The contribution from infinity vanishes for this shift, Inf = 0. This behavior for the full 
amplitude can be inferred from the corresponding behavior of the known one-loop QCD 
amplitude found by deleting cp, ^4 (1^, 2+, 3^, 4q). The rational part of this amplitude is 



a constant times the tree amplitude |Q, and it is easy to see from eq. ( |A.4| ) that the tree 
amplitude vanishes under the [4, 2) shift. Injecting a finite amount of momentum through 
the field (p should not affect the large- 2: behavior. We confirm this assumption a posteriori 
by checking factorization limits that are independent of the ones used to construct the 
recursion relation. It is also easy to verify that C4{z) vanishes at infinity; i.e., Inf C4 = 0. 

Next we look at the recursive diagrams R^. The only diagrams that give non- vanishing 



contributions are the ones shown in figure 10. We evaluate each of them separately. 



Diagram (a) is given by. 



D 



(a) 



^3(1 



2+ 



S12 



3 ('^>3t,4, 



(3.45) 



The tree amplitude yl3'^^(i;^, 3^, 4q, P^) is a simple MHV 0-amplitude. The loop three- 
point rational part R^^l^ , 2^, —P~^) can be extracted from the rational part of a one-loop 
splitting amplitude for g —>■ qq. It is equal to the MHV tree amplitude A^^\l'^ , 2+, —P^), 
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multiplied by the {0{e^)) rational part of of the loop splitting factor r^^ (it, q^,q^) defined 
in ref. |7l, 



(3.46) 



rat. 18 6 

Here 5r is a regularization-scheme dependent parameter, which fixes the number of helicity 
states of the gluons running in the loop to (4 — For the 't Hooft-Veltman scheme |87| 



Sr = 1, while in the four-dimensional helicity (FDH) scheme [S8, 8£] = 0. 
Thus we get for diagram (a), 



D 



(a) 



,[2(-P)]2\ /83 5_ 



[12] 



18 6 / su 



.(P4)2 
' (3 4) 



To remove the dependence on P we use the on-shell condition, 



(12) = ^ (12)+z(14)=0 ^ z 



(12) 
(14) 



and 



P = |1](1| + |2](2| = |l](l| + |2](2|+z|2](4] = P + z|2](4], 
plus some simple spinor product algebra, to get 



D 



(12) (3 4) V 18 6 J 4 Q' \ -.o 



18 6 



(3.47) 



(3.48) 
(3.49) 

(3.50) 



For diagram (b) we have the same on-shell condition as for (a). We also need the 
rational part of the leading-color (pqqg amplitude A^\4i, 1^, 2+, 3^). This can be extracted 
from the Hqqg amplitude |48| and the finite amplitude ^3^^((/), 1^, 2+, 3+) ||5^. Then a 
calculation very similar to that for diagram (a) yields 



D 



(b) 



A'^^ (l,r , 2+ -P+ ) — i?3 (0, 3 J , 4q , P- ) 



' [12] 



Sl2 



,(P4)2 
' (34) 



83 5n 

2 + 

18 6 



A 



(0), 



!.,lr,2+,3t,4-) X I — - 



(3.51) 



For diagram (c), the required one- loop (j) amplitude is A^^^ (0, 1^ , 2+, 3+) [|5§. The 



on-shell condition becomes 



[34] = ^ [34]-z[32]=0 ^ z 



[34] 
[3 2] 



and 



P = |3](3| + |4](4| = |3](3| + |4](4|-z|2](4] = P-z|2](4] 



(3.52) 
(3.53) 
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The diagram evaluates to 



By substituting z with its value given by the on-shell condition for P, eq. ( 3.52 ), and using 
the Schouten identity, we find that 

^(e) ^ , [2 3]^ / ^ ^ 1 [12](1|(2 + 4)|3] \ 



[12] [34] V 2 [2 3]S234 y 

- A^\6^ 1- 2+ 3+ 4-) ( 2 I 1 [12](1|(2 + 4)|3] 
,l,-,2^,3^,4Qj ^ ^2^^^^^^ 

In principle, there could be recursive diagrams associated with the S123 and S341 chan- 
nels. However, these diagrams vanish because on one side of the pole is a 4>qq amplitude. 
The amplitude A2{4': > 2|j') vanishes by angular momentum conservation, while the ampli- 
tude ^2(1^) ) 2^) vanishes because the quarks are massless and interact only with gluons, 
via chirality-preserving interactions. Similarly, the 523 and S41 poles are absent because 
there is no three-point amplitude containing two different flavor quarks. The sum of the 
recursive diagrams is 

< = Z)f)+I)f) + Z)f\ (3.56) 
Next we evaluate the overlap terms O4. They are given by 

0,= y Res^^^^, (3.57) 

poles a 

where the sum is only over the physical poles. In our case, physical poles can arise only when 
the following intermediate momenta go on shell (the same channels that admit possible 
recursive diagrams): 

(12) (2 3) , , 

p?, = o ^ ^ = -|T4y' ^23 = ^ ^ = ^' (^-^^^ 
Pi4 = o ^ ^ = p|' ^41 = ^ ^ = |rf' (3-^9) 

^^•^ (4|(1 + 3)|2]' (4|(1 + 3)|2] ^ ^ 

However, we note that the S23, S41, S123 and S341 channels had no recursive diagrams. This 
does not necessarily imply the absence of an overlap diagram (in principle CR/^ could have 
a worse behavior than i?4 in a given channel), but it is easy to check that eq. (3.42) has no 
poles in these channels. 

The two remaining cases correspond schematically to diagrams (a) and (b) shown in 
figure 0. We write 

O4 = O'l^ + of ^ , (3.61) 
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2+ 




2+ 




3+ 



(a) 



Figure 11: Schematic diagrams corresponding to the overlap terms O4 of A^^{<j),lg ,2+,3^,4q). 

The diagrams are computed from the residues of CRi{z)/z at poles in the z plane satisfying (a) 
P^^ = and (b) PI4 = 0. 



where diagram (a) corresponds to = and diagram (b) to P^^ = 0. It is simple to see 



from the [4, 2) shift of Ci?4 in eq. (3.42) that O^^^ also vanishes, because there is no factor 



of (12) in the denominator of Ci?4. Thus the only non- vanishing overlap contribution 
comes from 0^^\ due to the factor of [34] in the denominator of CR4^. The residue is easily 
found to be 



Oi = 



(b) ^ i (l|(2 + 4)|3] [23] 

2 [3 4]S234 



(3.62) 



We are finally ready to assemble the remaining rational terms of ^4 ((/), Ig , 2+, 3^, 4q), 
given by the sum 

^ ' ^ (3.63) 



i?4 — i?4 + O4 . 



Using eqs. ( |3.56| ) and ( 3.62 ), we get 



Ra = A 



(0), 



l^-%V24^0^1.-,2^3t,45), 



9 



3 J 



(3.64) 



which is also recorded in eq. (4^). Note that the term with the unphysical pole in S234 in 
the recursive diagram D^-^ cancels against the overlap diagram of ^ for the same channel. 
Such recursive-overlap cancellations are a common feature. 



3.4.2 Spurious poles in rational parts 



The cut-completion process described in eqs. ( p. 34 ) and ( 3.35| ) removes certain types of 
spurious singularities, namely those associated with terms like (lnr)/(l — r)" for n > 1, 
where r is a ratio of two momentum invariants that differ by one massless external leg. It 
ensures that C„(z), and therefore Rn{z), is free of spurious singularities when r{z) — > 1, 



where r{z) is the shifted value of r. In refs. [56, 57, 51[ it was found that this type of 



cut completion was sufficient to remove all spurious poles in the z plane for a large class of 
QCD and (p amplitudes. (Note that if a spurious denominator factor is unaffected by the 
particular shift used, eq. ( p. 271) , then it will not produce a spurious pole in the z plane.) 
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However, poles of the type r{z) = 1 certainly do not exhaust the set of potential spurious 
poles for a general amplitude. These poles correspond specifically to Gram determinants 
associated with two-mass triangle integrals. In a general amplitude, there are also poles 
associated with the Gram determinants of a multitude of scalar box and triangle integrals 
with varying numbers of external masses. It would be very difficult to construct a cut 
completion C„ that removed all spurious poles. 

On the other hand, a completely general alternative method for handling the spurious 
poles was sketched in ref. ||5^, and fully implemented numerically in the BlackHat pro- 



gram |44|. This method did not use cut completion, but rather the original decomposition 



An^ = Cn + Rn, and instead relied on the fact that the residues of Cn{z)/ z and Rn{z)/z 
cancel at every spurious pole. The contour integral of Rn{z)/z at infinity requires the sum 
over spurious pole residues of Rn{z)/z, but this sum can be evaluated using the shifted cut 
part Cn{z), as 

_ y Res^= y Res^. (3.65) 

spurious poles 13 spurious poles /3 

Here we simply note that a hybrid approach is also feasible: First one removes the 
spurious poles that can be easily removed by cut completion, such as eqs. ( |331 ) and (|33^ ). 



This procedure leads to the overlap terms ( |3.4[1|) in the usual way. Then one considers the 
contour integral at infinity of the remaining rational terms Rn{z) / z. One evaluates the sum 
over residues of Rn{z)/z, at the spurious poles z^ that were not removed by cut completion, 
by using the fact that their residues still cancel against those of the completed cut terms, 
Cn{z)/z: 

_ y Res^= y Res^, (3.66) 

spurious poles 7 spurious poles 7 

where {7} is a subset of {(3}. 

Note that Cn{z) in eq. ( |3.66| ) includes rational terms as well as cut terms; whereas 
Cn{z) in eq. ( |3.65| ) is a pure cut function, which nevertheless can have rational- function 
spurious-pole residues after Taylor expansion around the pole. In fact, it is only the 
rational- function part of the residue of Cn{z)/z that we require; the terms containing 
logarithms, polylogarithms and vr^ factors are guaranteed to cancel, because the residue of 
the rational function Rn{z)/z can have no such terms. 

In our calculation of the (pqqgg amplitudes, after removing the r{z) = 1 spurious poles, 
we found that certain spurious poles still remained, due to denominator factors in integral 
coefficients of the form (ij) or [ij], where i and j are color non-adjacent legs. If i and j 
are color-adjacent, then (ij) and [ij] denominator factors represent physical singularities, 
corresponding (for real momenta) to the region where ki and kj become collinear, Sij — > 0. 
However, physical collinear poles for color-ordered primitive amplitudes only occur when 
i and j are color- adjacent; in the non-adjacent case, the {ij) and [ij] factors generate 
spurious singularities. They are associated with the "easy two mass" box integral in which 
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the two diagonally opposite massless momenta are ki and kj, whose Gram determinant 
contains a factor of Sij = (ij) [ji]. The factors (ij) and [ij] can also be seen in the 
denominators of the solutions for the loop momenta in the quadruple cut for the easy-two- 
mass box kinematics (see e.g. eq. (2.7) of ref. p^]). The same factors also persist in the 
limiting case of a one-mass box. 

Factors of {ij) or [ij] do appear in the denominators of coefficients of scalar integrals 
in the QCD and (j) MHV n-gluon amplitudes with non-adjacent negative-helicity gluons 
labeled i and j, whose rational parts were computed using on-shell recursion relations |58 
52]. However, in these cases the [i,j) shift was used, which leaves both (ij) and [ij 
unshifted, and therefore produces no spurious pole in this channel. 



In this work we encountered spurious poles associated with {ij) factors in several 
amplitudes. Here we illustrate how to use the method described above for the specific 
example of the leading-color primitive amplitude ^4 (</>, 1^, 2+ 3+, 4~). Using the methods 



of sections 3.1-^^, we have calculated the cut-containing parts of ^4^ (</>, 1^, 2+, 3+, 4~), 
and obtained the completed-cut terms C4, which can be read off eq. ( [4. 191 ) by ignoring the 
purely rational terms in the last few lines. To compute the rational part of this amplitude 
we chose to use a [4, 1) shift. Inspecting 6*4, we see that several terms contain the spurious 
denominator factor (13), which will potentially lead to a spurious pole under the [4,1) 
shift. These terms are. 



(14)^ 



(12) (34) (13) 
i (12)2 [2 3]^ (3 4) 



Ls_i (si2, S23; S123) + Ls_i (S34, S41; S341) 



^123 



(13) 



^^_{_^2) (34) [23]2(14) Li 



'12 



(13) 



'12 



(3.67) 



The spurious pole satisfies = (1 3) = (1 3) + 2; (4 3), or z = (1 3) / (3 4) = Zsp 
Using eq. ( |3.66|) , we need to compute 



Res 



(3.68) 

This is because the 

relevant scalar box integral m D = 6 dimensions, which is nonsingular as S13 — > 0, can be 
written as 



Now the Ls_i functions in eq. ( 3.67 ) actually vanish as z 



rD= 



\S12, S23; S123) 



-I cr ■ 



LS- i(gi2,g23;'Sl23) 
•S13 



(3.69) 



and similarly for the other Ls_i function. Because T^^^ is smooth in this limit, the Ls_i 
functions must contain a factor of S13 = (13) [31] in the limit S13 0. Thus the terms 
containing the Ls_i functions in eq. ( 3.671) do not contribute to the residue. 



After expanding the remaining logarithms and rational terms in eq. ( |3.67| ) around 
z = Zsp, we find that the logarithmic part of the residue cancels, as expected. Keeping the 
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rational part of the residue, and simplifying, we get 



C,{z) _ i (34) (14)^ [23] .2 (14) [23] (34) 

^=lp z 6(2 3)(13)(4|(1 + 3)|2] ^ 3 [1 2] (2 3) (1 3) ' ^ ' 

This term has to be added to the recursive diagrams, overlap terms and CR^ to complete 
the full rational terms of A^{(f), 1^, 2^, 3"'', 4~). The full rational terms, as well as the full 
amplitude, are now free of spurious singularities as (13) — > 0. 

The procedure outlined above can be performed in a systematic way for any amplitude, 
since the locations of the possible spurious poles under a chosen shift are known a priori, 
or they can be inferred simply by inspecting the completed-cut terms, Cn- Whenever, after 
absorbing spurious singularities according to eqs. ( 13.34] ) and ( p.35 ), we are left with residual 



spurious poles, we can always compute their contribution to the remaining rational terms 
by evaluating the corresponding residues of Cn{z)/z instead. 



4. The one-loop HqqQQ and Hqqgg amplitudes 

In this section we present our main results for the one-loop HqqQQ and Hqqg^g^ ampli- 
tudes. First we outline how to obtain all primitive ^-amplitudes, using only a minimum 
set of them. Then we give the full analytic expressions for these amplitudes, followed by 
numerical results at a specific kinematic point. We then show how to obtain the color- and 
helicity-summed cross section for a pseudoscalar Higgs plus two quarks and two gluons. 



using our results and those of ref. [21]. As another application, we show how to compute 
part of the virtual one-loop color-singlet interference term between the gluon-fusion and 
VBF Higgs production mechanisms. Finally, we mention the various consistency checks we 
used to verify the correctness of our expressions. 



4.1 Preliminaries 

We obtain the one-loop corrections to Ai{H, 1^, 2g, 3q, 4q) and Aa{H, Ig, 2g, 3, 4) by com- 
puting color-ordered primitive ^-amplitudes in a helicity basis, following our discussion in 



section 2.2. Once we have the complete set of ^-amplitudes, the (/)T-amplitudes are obtained 



by parity, eq. (|2? 

Consider first ^4((/), 1^, 2q, 3q, 4q). Because the quarks are massless, chirality is pre- 
served along a quark line. By convention, all external legs are outgoing, so the helicities of 
any quark-antiquark pair have to be opposite. Thus we need only consider the four helicity 
configurations ^4((/>, 1^"^, 2^ , 3^^, 4q), where A,A = it are the helicities of the q and Q 
quarks, respectively. 

Suppose the anti-quark q (leg 1) has positive helicity. We can obtain these cases from 
the cases where it has negative helicity by using charge conjugation, which reverses both 
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quark lines {q ^ q and Q ^ Q): 



^4(</', 1^, 2-, 3^, 4+) = A4{4>, 2r, 1+ 4t 3^ 



^4(</', 2^,1+ 4^,3+) 



(4.1) 
(4.2) 



^4(<?^,1^2-,3t4Q 

Now taking the anti-quark q to have negative hehcity, we see that there are two independent 
hehcity configurations that we need to compute, 



A^i^, ir,2+ 3t 4^ 



and 



Aiicj), ir,2+ 3;=,4+). 



(4.3) 



Here A4 is shorthand for the three types of primitive amphtude in this case (Ic, sic, and f). 
Recall from eq. (p.9|) that parity gives the i;^>t-amplitudes in terms of the (/^-amplitudes. 



^4(</'^l^2-^3J,4, 



Aa 



(4.4) 



where the operation (ij) <^ [j i] conjugates spinors but does not reverse the sign of ab- 
sorptive parts of loop integrals. 

For the Hqqgg amplitude, there are two cases to consider, depending on whether the 
helicities of the two gluons are the same or opposite. In the case that they are the same. 



say both positive, we have, using the decomposition (2.7) and parity, 
^4 (^^', 1,- , 2+ , 3+ , 4+ ) = ^4 (0, 15 , 2+ , 3+ , 4+ ) - U4 (0, 1^ 2- , 3- , 4- ) 



(4.5) 



The amplitude ^4((/), 1^ , 2+, S"*", 4"*") vanishes at tree level. For this reason, the one-loop 
amplitude is quite simple and is given below in eqs. ( ^l6|) , ( FTtI) and (p^ ). However, 
the amplitude ^4(0, 1^, 2~, 3~, 4~) is next-to-maximally-helicity violating (NMHV), and at 
one-loop it is considerably more complex. (For example, the coefficients of the three-mass 
triangle integrals are nonzero for this amplitude.) We will leave its analytic computation 
for future work. 

Instead we turn to the case of opposite-helicity gluons, which can be decomposed as 
^4(ii',l5,2+,3±,4^) =^4(<^,l5,2+,3±,4T)- U4(</', 1^ 2^^, 3^, 4±)] . (4.6) 

In this case the (j) amplitudes are both MHV, and of a similar complexity as the (pqqQQ 
amplitudes. Again, using charge conjugation we can exchange the roles of anti-quark and 
quark, so as to obtain the remaining i?i-amplitude helicity configurations, in which the 
anti-quark q has positive helicity, 

^4(1^ 2-, 3±, 4^) = ^4(2,-, 1+ 4^, 3±). (4.7) 

Using the color decompositions in section |2.2| , the problem is reduced to computing the 
primitive amplitudes ^4 , Af^ and ^4 for two four-quark helicity configurations, 1^2^3^4^, 
and A^, Af and A\ for two two-quark- two-gluon helicity configurations, and two color 
orderings, namely 1^2+3^4^ and 1^3^2+4^. 
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4.2 Full results 



In section |^ we showed in specific examples how to compute various ingredients necessary 
to obtain the full ^4 1^, 2+, 3^, 4g) and A^{4>, 1^, 2+, 3+, 4") amplitudes. We used the 
same techniques for the quadruple cuts and ordinary two-particle cuts in all channels, and 
for all other color components and independent helicity configurations, in order to arrive 
at the full results for the ^-amplitudes. 

It is worth noting that the computation of the (j)qqQQ primitive amplitudes in both 
helicity configurations was significantly simpler than that of the (pqqgg amplitudes. The 
expressions were more compact at each stage (due in part to the higher symmetry of these 
amplitudes). In addition, we had no remaining spurious poles for the shifts we chose 
to perform. Specifically, we used a [4,2) shift for all the ^4((;/), 1^, 2+, 3^, 4g) primitive 

amplitudes, and encountered no contributions from z 00 (i.e., 6*4(2) as well as A^\z) 
vanish in this limit). We used a [1,3) shift for A4 (</), 1^, 2+, 3^, 4q). Here there was 

a contribution from 2; — > 00, from C4. We used a symmetry to obtain the sic and f 
components of this helicity configuration from the previous one. 

The (j)qqgg amplitudes were more intricate. For A4((/), 1^, 2^, 3"'", 4~), we used the [4, 1) 
shift, and we had to compute residues of spurious poles in the L, R and fermion loop (f ) am- 
plitudes, although there were no contributions from 2; ^ 00. For j44((/), 1^, 2+, 3~, 4+), we 
used the [3, 2) shift. There were not only spurious pole residues in the L and R components, 
but also a contribution from z ^ co (from C4) in the L component. We cross-checked our 
results for the L and R components of ^4((/>, 1^, 2^, 3^, 4^) using the [2,4) shift. 

Finally, for the (j)qgqg L amplitudes, we used a [4,1) shift for ^^((/>, 1^, 2+, 3+, 4~), 
and a [2,1) shift for (0, 1^ , 2", 3+, 4+). There were neither spurious pole contributions, 
nor contributions from z ^ 00. The full results, after assembly and simplification, are 
presented below. 



4.2.1 cj)qqQQ 



For the (pqqQQ ( — \--\ — ) configuration we have 



-i4^(0, i5,2+,3t,4-; 

-i(12)(34) [23]^'^' 



+ 2i4°)(</)T,lr,2+,3t,4Q), 



-iA 



(0), 



l5'2+,3t,4Q)xy" 



-^123 
-Sl2 



Li 



+ 



-5234 
-S34 



'12 



'34 



2(14) [2 3] 



-S12 



S12 



Lo(- 
+ ^ 



-5234 
-534 



S34 



(4.8) 
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with 




- Ls^J^^ (si23, S234; S23, mff) - Ls^J^^ (5341,5412; S4i,mjj) 

- Ls_i (S23,S34;S234) - Ls_i (§34, S41 ; S34l) 

- Ls_i (s4i, S12; S412) - Ls_i (si2, S23; S123) + ^ - , (4.9) 



-af(</>,ir,2+ 3t,4Q) 



-iA 



(0), 



ir,2+ 3+ 45)xF 



sic 



+ -(12) (3 4) [2 3]^ 



Li (^)^Li (^) 



'12 



*34 



(14) [2 3] 



S12 



S34 



(4.10) 



with 



V 



sic 



2 \ ^ 



-S12 



+ 



2 \ en 



-S34 



_3_ 
2^ 



2 \ e 



-S12 



+ 



2 \ e 



-S34 



- Lsi^'' (s412,Sl23;Sl2,m^f) - LsJ^'^ (§234, S34I ; ^34, - 7 - Sr , (4.11) 



and 



4 (,^, ir , 2+ 3 J , 4- ) = (^, ir , 2+ 3t , 4- , . - - 



/V, 



2 \ 



-S12 



+ 



,2 \ 



-S34 



20 



(4.12) 



For the (pqqQQ ( — I h) case, we find 

-z4(</.,ir,2+ 3^,4+) = -af (</.,l,-,2+ 3^,4+) 

2" 



1 (12) (3 4) 
2 (24)2 



/(14H23)Y 
V(13) (24); 



Ls_i (S23, 534; S234) + LS_1 (S4i, S12; S412) 



Jl 



'12 



'34 



-S34 



+ 2 



[24] 
(24) 



(23) 
[14] 



T / — £412 1 
S2A ^ 



S12 



+ ln 



( -S412 \ 
V -^12 J 



+ 



(14) 
[2 3] 



S24' 



"0 I, -534 ) 



S34 



+ ln 



S34 



+ 



(14) (2 3) [2 4] 



(2 4) 



Lo(^^Lo(-^ 



S41 



S23 



1 1 
2(24p 



(34)(S24-S4l) , (12)(S24-S23) 



[12] 



+ 



[34] 



+ 2z4°)(<^t^l-^2+ 3^,4+). (4.13) 
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The sic and f primitive amplitudes for the (j)qqQQ (— H h) are simply related to those for 

(— +H — ), because two of the external legs can be exchanged at the cost of a minus sign, 



A 



sic / 



^ 1 2^' ^Q' '^o) 



4(<^, ir,2+,3^,4+) 



f( 

i 1^ 



(4.14) 
(4.15) 



These relations were used already in constructing the partial amplitudes ( 2.14| )-( 2.17] ). 



The result for any color or helicity component of A^I\h, Ig, 2^, 3q, 4q) can be readily 
obtained using eqs. (^j])- 



4.2.2 cj)qqgg 



The results for the infrared- and ultraviolet-finite helicity amplitude 4>qqgg (—+++) can 



be extracted from ref. |5C]. We give them here for completeness: 



,lr,2+ 3+,4+) 



1 (12)(1|(3 + 4)|2] 1 (13) [34] 
2 (2 3) (34) (41) 2 (2 3) (34) 



+ 2 



(1|(3 + 4)|2]^ 



(1|(2 + 3)|4]2(2|(1 + 3)|4] 



(34)(41)(3|(1 + 4)|2] (12)(23) (3|(1 + 2)|4] 



[2 4]- 



m 



H 



1(13) [3 4] (41) 



[12] S4i2(3|(l + 2)|4](3|(1+4)|2] 3 (12)(34)' 



(4.16) 



-af(,^,ir,2+ 3+,4+) 



(1|(2 + 3)|4] , (12) [2 3] (31) 



+ 



(23) (34) (23) (34) (41)_ 



(4.17) 



-iA{{ 



ir,2+,3+,4+) 



1 (13) [34] (41) 
3 (12) (34)^ 



(4.18) 
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For (f)qqgg ( — \—\ — ) we obtain, 



-ai'(0,ir,2+ 3+,4-) = -af (</>,ir,2+ 3+,4 



- Ls_i (s23, •S34; S234) - Ls_i (s4i, S12; S412) 



+ 



+ 



+ 



(14^ 



(12) (3 4) (13) 
4(13)2(4|(1 + 2)|3]3 



Ls_i (si2, S23; S123) + Ls_i (S34, S41; S341) 

1 (12)2 j23]^ (34)2 



(12) (34) 



-(12) [23]^ (34) (4|(1+2)|3]- 



(13) 



-S123 

-Sl2 . 



1 (13)2(24)(4|(1+2)|3]2 (13)(14)(4|(1+2)|3]2 1 (1 2)(34)[2 3]2(14) 
2 (12) (23) (34) (12) (34) 2 (13) 



*12 



'12 



l(12)(34)(2 4)[2 3]2Ln^j , (1 4)^ (4|(1 + 2)|3] 



Ln 



2 (2 3) 

0(1 4) (2 4) [2 3] 



+ 



-^123 



'34 



(12) (34) 



S12 



Ln 



- 



A' 123 



S12 



Lr 



Sl2 



+ 



-^234 
-S34 



S34 



(14)^ 



+ 



(2 3) 

-2iAf\(t>, ir,2+ 3+,4- 

5 (14)2(4|(1 + 2)|3] _ 1 (14)2 [2 3] (3 4) 

6 si2(12)(34) 



(12) (3 4) (13) 



In 



-S123 
-S12 



2 (14) [2 3] (34) 
6(23)(13)(4|(1 + 3)|2] ' 3 [1 2] (2 3) (13) 



+ 



2(14)(24)(4|(1 + 3)|2] , 1 [23] (4|(1 + 3)|2](4|(2 + 3)|1] 
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4.2.3 (}>qgqg 

For (pqgqg, we have the "reflection" relation ^4^(lg, 3, 2q, 4) = A^{lq, 4, 2q, 3), so we do not 
need to quote Af separately. 

Again we take the results for the infrared- and ultraviolet-finite helicity amplitude 
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The results for (pqgqg ( — \—\ — ) are given by 
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The results for (pqgqg ( are 
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Using eqs. (|2.g| ) and (4/7), one can obtain any color or helicity component of the 



A^I\h, lg,2q,3=^,4^) amplitudes. 



4.3 Numerical results 



In order to facilitate comparisons with future work, we present here numerical values, at 
a single phase-space point, for the bare, unrenormalized primitive amplitudes computed 



in the paper. We choose the same kinematic point as in ref. |21|, the configuration H 
Q1Q2Q3Q4 in which the Higgs (or (p) and parton four-momenta take the values 



= (-1.0000000000, 0.00000000000, 0.00000000000, 0.00000000000), 
ki = (0.30674037867, -0.17738694693, -0.01664472021, -0.24969277974), 
k2 = (0.34445032281, 0.14635282800,-0.10707762397, 0.29285022975), 
A:3 = (0.22091667641, 0.08911915938, 0.19733901856, 0.04380941793), 
fc4 = (0.12789262211, -0.05808504045, -0.07361667438, -0.08696686795). 



(4.35) 



We substitute // = mn = 1, and use the 't Hooft-Veltman scheme [37|, in which 6r = 1 
(in accord with ref. ||2l|). Discussions of the conversion between different dimensional 



regularization schemes can be found in refs. |8G, 89, 21]. The dependence on in our 
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Table 1: Numerical values of cfyqqQQ primitive amplitudes at kinematic point (4.35) 
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Table 2: Numerical values of ipqqgg primitive amplitudes at kinematic point (4.35) 



formulae agrees with the shift in HqqQQ and Hqqgg amplitudes between HV and FDH 
regularization schemes that is quoted in ref. [pl| ]. 

In tables ||, ^ and ^ we present numerical values for the unrenormalized primitive 
amplitudes computed in the paper. Note that the overall phases are convention-dependent. 
Phase-independent quantities can be constructed by dividing by the corresponding tree 
amplitude. The tree amplitude is identified as —1/2 of the sic entry in table the 
negative of the R entry in table ||, and —1/2 of the L entry in table ^. 

In table Q we give the numerical value of the virtual correction to the color- and 
helicity-summed cross section for the HqqQQ process, according to eq. (|2.18|) but omitting 
an overall factor of 2C'^crg^{N^ — ^)Nc. The result is constructed from the one-loop 
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Table 3: Numerical values of (pqgqg primitive amplitudes at kinematic point (4.35) 
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Table 4: Numerical value of the one-loop correction to the HqqQQ cross section at kinematic 
point (|4.35|), omitting an overaU factor of 2C'^crg^{Nl ~ l)Nc from eq. (|2.18|). 



amplitudes given in this paper and the tree amplitudes ( |A.4| ). The dependence on the 
number of colors Nc and massless quark flavors n./ is shown explicitly. If one substitutes 
Nc = 3 and Uf = 5, adds the contributions, and multiplies by 1/4 x (A'^^ — l)Nc, then the 



result agrees with that for process A in table I of ref. |21]. (The factor of 1/4 arises because 
a factor of = (2C)^ is extracted instead of in ref. 121 



To convert the bare, unrenormalized amplitudes presented here to renormalized ones 
in an MS-type subtraction scheme, one should subtract the quantity 
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(4.36) 



from the corresponding one-loop amplitude A\ in eq. dmi ) or ( ^^201 ). After this subt rac- 



tion, the rational parts of the 1/e poles are purely infrared, and take the form of a sum 
over contributions from each external parton. 
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(4vr) 
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1 /II 
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for (pqqQQ, 
for (l)qqgg, 



(4.37) 
(4.38) 



in accordance with the general form of infrared singularities of one-loop amplitudes |90]. 



A final contribution that needs to be included is the one-loop correction to the Hgg 
effective operator of eq. (p.lj) , which shifts its coefficient from C to C x [1 + llas/(47r)] p3|| . 
At the level of the NLO virtual cross section, this effect can be taken into account by an 
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addition to eqs. ( 2.11 ) and ( 2.2C| ) of the form 



11 j^,aT . (4.39) 

4.4 Pseudoscalar Higgs amplitudes and cross section for Aqqgg 

As a byproduct of our calculation of the (/>-amplitudes, we can obtain the respective am- 
phtudes where the scalar Higgs boson H has been replaced by a pseudoscalar Higgs boson, 
A. Pseudoscalar Higgs bosons are present in many extensions of the SM, such as the 
MSSM. Here we assume that we are in a kinematic regime where the production of the 
A boson plus jets can also be treated in the large mt limit. As mentioned earlier, the 
overall constant C is different for the A case [^0[. Otherwise, the only difference between 
the two computations is that instead of taking the sum of the (/>- and (f)^- components, the 
pseudoscalar amplitudes are given by their difference divided by i, according to eq. (|2.8p. 



Furthermore, we can combine our results with those of Ellis, Giele and Zanderighi 
(EGZ) [^l|, to obtain the color- and helicity-summed cross section for Aqqgg. EGZ used a 
semi-numerical method to compute the color- and helicity-summed cross section for Hqqgg, 
which we may write schematically as, 

cTEGZ = E [^2)*(A)^«(A) +^g)(A)^«*(A)] = 2Re [-4g^*(A)^«(A)] | . (4.40) 

Here by AhW we denote a Hqqgg amplitude in a helicity configuration A, and the sum 
is over all possible helicity configurations. In every summation in this subsection there 
is also an implicit sum over colors which is given by eq. ( |2.23| ), but for simplicity we do 
not write it out here. The helicity sum includes the two independent MHV configurations 

( — hH — ) and ( — I h) analyzed in this paper, and their conjugate ones (obtained by 

parity). Moreover, they include configurations ( — I ), ( — 1-++) and their conjugates. 

These last configurations require NMHV (j) amplitudes, which we did not compute here, 
and which would be needed to provide a full analytic description of the one-loop Hqqgg 
and Aqqgg helicity amplitudes. However, we will see that because the contribution of the 
latter configurations to the Hqqgg cross section is encoded in ctegz, it is possible to use 
this result instead for the purpose of computing the Aqqgg cross section. 



Note that we can rewrite eq. (|4.4C| ) as 

CTEGZ = 4Re ( V [Ag)* (AO A« (A')! | , (4.41) 



with A' now labelling each of the four helicity configurations with fixed helicities for the 
antiquark-quark pair, q~q'^, namely ( — 1-±±). Also, from eq. ( [4.41 ) we see that cjegz 
contains the NMHV Hqqgg amplitudes in the quantity 

a™^ = 4Re{A2^*(A_)Ag)(A_)+^S)*(A+)^«(A+)} C ctegz, (4.42) 
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with A_ = ( — I ) and A+ = ( — We can extract ct^q^^ from eq. ( 4.41 ), know- 
ing the Hqqgg MHV amphtudes inferred from our formula (4.16)-(4.34). Hence we will 
consider <7egz^^ known from now on, and we will use it in order to compute the color- 
and helicity-summed cross section for Aqqgg, a a- Similarly to eqs. ( 4.40 ) and ( [4.41| ), a a 
is given by 



<TA = 2 Re 




■(A)4^(A) 



4 Re 




(4.43) 



with Aa{X) denoting the pseudoscalar Aqqgg amplitudes. 



Because a a is expressed as a sum over A', we focus on the four configurations ( — hH — ), 

( — I h) and ( — I ), ( — We can straightforwardly compute the terms coming 

from the first two (MHV) configurations using our results from section |4.2| and eq. (^. j), 
as mentioned in the beginning of this section. For the last two (NMHV) configurations, we 
need to compute the quantity 
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Let's look at each amplitude in this expression separately. The tree amplitudes are simple. 



because A 



and (|2.8D, we have 



(0) _ (^(0) _ ^W)/^, and ^J^t^(A-) = A^^\x+) = 0. Therefore, from eqs. (^) 



4^A-) = i<)(A_) = i4^(A. 



<(A+) 



1 



-Ag)(A+), 



(4.45) 
(4.46) 



We choose to express the one-loop amplitudes A^-'(A_) and A^'(A+) using eqs. 
and ( ^.81 ) in the following way 



4^(A.) = yK^(A.)-A^V(A.)^ 



^iJ)(A-)-2^W(A- 



A«(A+) + 24^^(A+) 



(4.47) 
(4.48) 



Substituting eqs. (|1|)-(|^) into (|^) we find that 



_NMHV 



4Re< 



>S'(A-: 



I L 



A«(A_)-2A^V(A-) 



+ 



-4Ha^) 



1 r 



^«(A+) + 24^)(A+ 



4Re{4)*(A_)4)(A_) + <)*(A+)4)(A+)} 

-8Re{4)*(A_)A(V(A_)+^2)*(A+)A«(A+)} 
«^ - 8Re {<*(A_)A«(A_) + (A+)<)(A+)} 



(4.49) 
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We notice that the first term in eq. ( 4.4S| ) is given by eq. ( [4.42| ), and the second term 
contains only tree and finite one-loop helicity amplitudes with cf) and (f)^ . The required 
NMHV tree amplitudes (see e.g. refs. ||, |^]) are given in eqs. ( [Al6| ), ([A?7|), ( |A^ and 
( |A.1C| ). The finite one-loop amplitudes are given in eqs. ( |4.16D -( |4T8| ) and ( [4.27| )"( p^ ). 
Therefore, we know all the ingredients necessary to obtain the full color- and helicity- 



summed cross section for Aqqgg, albeit only semi-numerically. Our results of section 4.2 
can be used to convert cjegz into aA- 



4.5 Interference with VBF production 



Our amplitudes for HqqQQ can be used to calculate analytically part of the interference 
between the qQ HqQ gluon fusion process and the tree-level vector boson fusion pro- 
cesses. Both these processes have the same initial and final states. However, at tree level 
one has a color-octet exchange and the other a color-singlet exchange. So there is no inter- 
ference at tree level. (For identical quarks, the exchange term does produce an interference, 
but it is extremely small [|60|].) At one loop, however, the color-singlet part of the one-loop 
correction to the gluon-fusion HqqQQ amplitude can interfere with the tree-level VBF 
amplitude, and we will provide an analytic formula for this contribution. This is only part 
of the virtual correction; the other part comes from the interference between the tree-level 
gluon fusion and one-loop VBF HqqQQ amplitudes. In addition, there is a real correction. 



The sum of all three terms has been computed numerically in refs. |24, 25| and it is quite 
small. 



We obtain the color-singlet part of HqqQQ from the corresponding (/>-amplitude, 
A4;s((/>, Ig, 2q, 3q, 4q), using eq. (2^). From the color decomposition ( p. 11 ), contracted 



with 5^^^ , we see that the color-singlet part is given by 



A4;s((/', Ig, 2g, 3qAq) = ^4;i(</', Ig, 2^, 3q, 4q) + A^,2{^, Ig, 2g, 3q, 4q). (4.50) 
Using eqs. (|2.14D -( ^rT7|) we get 



A. 



4;sl 



1- 9+ o+ A- 



- 1 



A!f(0,l5,2+,3t,4Q)+4^(0,ir,2+,4^,3+) , (4.51) 



^4;s(</',l,-,2+ 3^,4+) = ^^ [4^(</.,lr,2+,3^,4+)+4'=(</.,ir,2+,4+,3Q)J , (4.52) 



in terms of the leading-color primitive amplitude ^4 ((/>, 1^, 2g, 3q, 4q). The relevant tree- 
level HqqQQ VBF amplitudes involve only ZZ fusion, not WW or WZ] they are given 

by 



A 



(0) 

4;VBFV 



{H, Ig ,2+, 3^, 4^ 



2i- 



rrir 



(14) [2 3] 



V (si2 - m|)(s34 - m|)' 



^4;VBFV 



{^1 ^q ) 2^, 3q, 4g) — A^.yj^-p(H, Ig ,2^, 4^, 3, 



(0) 



(4.53) 
(4.54) 
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with V the vacuum expectation value of the Higgs field. Finally, the color-singlet virtual 
interference between the two processes is 



2alN^ Re 



"44°vBf(-^' ^9' ^9' 3q, 4Q)^4;s(i?, Iq, 2g, 3q, 4q) 



(4.55) 



Equation ( [4.55 ) is to be understood with an implicit summation over all allowed polar- 
ization states of the external quarks. We implemented eq. ( [4.55| ) numerically at a few 
phase-space points and obtained agreement |61] with this part of the full interference com- 



puted in ref. |24|. We did not perform such a comparison against ref. |25]. 



4.6 Consistency checks 

It is important to verify that the methods and the results presented in this paper yield 
the correct answers for the HqqQQ and Hqqgg amplitudes. We have used three types of 
independent and non-trivial checks on our expressions. They are based on collinear limits 
that the amplitudes should satisfy, symmetries under which they should remain invariant, 
and numerical comparisons with previously computed expressions. For Hqqgg only the 
first two types of checks were possible, whereas for HqqQQ all of them were performed, 
providing an even more solid check. We have found that our amplitudes agree with all the 
checks, and we outline the process further in the remainder of this section. 



4.6.1 Collinear behavior 

A powerful handle on the correctness of the amplitudes is their collinear behavior. When 
two neighboring external legs become collinear, an n-point amplitude has to correctly fac- 
torize onto an {n — l)-point amplitude, multiplied by the corresponding splitting amplitude 
for the two collinear legs. In the case of one-loop amplitudes, the factorization is onto a 
sum of possible factorizations with the loop belonging either to the splitting amplitude, or 
to the remaining {n — l)-point amplitude [^]. There is also a sum over the helicity h of the 
intermediate state P carrying momentum /cp ~ 0. For the (pqqQQ amplitudes, in the limit 
that momenta ki and k2 become collinear the factorization is onto a (pQQg amplitude, 



4^)(0,15,2„3q,4q) ^ Yl 



h=± 



4')(0,3q,4q,P^) X SplitL°i(l5,2,;z) 
+ 4°)(,/.,3q,4q,P'^) X SplitL'l(l,-,2,;z) 



(4.56) 



with kp = ki + k2, ki zkp, and k2 ~ {1 — z)kp. The splitting amplitudes depend 
on the longitudinal momentum fraction z, which is the momentum fraction carried by 
leg 1, a real variable with < z < 1. (It is unrelated to the complex variable z used 



for the shifts performed in the previous sections.) Replacing (j) with H in eq. ( 4.56 ), we 
get the collinear behavior of the Higgs amplitudes, while the 3 || 4 collinear limit can 
be obtained by exchanging q and Q. The 1 || 2 and 3 || 4 limits are the only collinear 
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limits of A^\(j), Iq, 2q, 3g, 4q) that exhibit universal singular behavior; there is no splitting 
amplitude for quarks of different flavor. 

The extension of eq. ( 4.56| ) to the 4'QQ99 amplitudes is straightforward. In this case, 
however, there are additional factorization channels, including channels where a gluon 
becomes collinear with an adjacent quark or gluon. For example, for the 2 || 3 and 3 || 4 
collinear limits we have 



A 



2 3 



E 

h=± 



A'i\(l,, l5,P,^4)xSplitL"i(2„3;z) 



+ <^(</.,l,-,Pg\4) X SplitL'^(2„3;z) 



',lg-,2„3,4) 



3 4 



E 

h=± 



4')(0, l5,2g,P^) xSplitL";:(3,4; 



(0), 



+ Af (,/.,l,-,2„P'^)xSplit^^^(3,4;z 



(1), 



(4.57) 



(4.58) 



and similarly for the Hqgqg amplitudes. Eqs. ( 4.56| )- (|4.5^ ) apply separately to the primi- 
tive amplitude components Ic, sic, L, R and f, after extracting the respective pieces of the 



one-loop three-parton amplitudes and splitting amplitudes [73|. 

We have checked that our expressions factorize correctly according to eqs. ([4.56|) - (|4.58|) 
and their analogues, for all possible non-trivial (singular) collinear limits. 



4.6.2 Symmetries 

The one-loop amplitudes we computed in this paper are MHV four-point amplitudes, with 
an equal number of positive- and negative-helicity external legs. As a consequence, they 
have to satisfy certain non-trivial symmetries, that become manifest only after assembling 
together all the pieces into a full answer. 

The symmetries of the four-quark HqqQQ primitive amplitudes can be summarized in 
the following form: 

• reflection or quark-exchange symmetry 

(i/, ir , 2+ , 3+ 4- ) = {H, 4r , 3+ , 2t , 1- ) , (4.59) 
(//, ir , 2+ , 3^ , 4+ ) = {H, 3r , 4+ , 1^ , 2+ ) , (4.60) 

• parity conjugation symmetry 

a['\h,1^,2^,3±Aq) = 4\h,2^,1^,4^,^q)\^.^^^^^.^, (4.61) 
4^)(i/,ir,2+ 3^,4+) = 4^)(i/,2r,l+,4^,3+)|^^^.^^^^.^^. (4.62) 
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The reflection symmetry properties ( [4.59D ~( [l.60 ) are satisfied by the component (/> and 0^ 
amphtudes as well, whereas the conjugation symmetry ( [4.61| )-( 4.62D only holds for the H 
amplitudes. 

For the two-quark-two-gluon Hqqgg and Hqgqg primitive MHV amplitudes, although 
there is no reflection symmetry, the following parity conjugation symmetries hold: 



Hqqgg conjugation symmetry 

A^I\h, ir,2+ 3+,4-) 
A^I\h, ir,2+,3-,4+) 

Hqgqg conjugation symmetry 

A^I\h, ir,2+,3+,4-) 
A^I\h, ir,2-,3+,4+) 



-4^)(//,2r,l+,4+,3-; 
-4')(i/,2r,l+,4-,3+: 



-4')(if,3r,4+,l+,2 
-4')(if,3,-,4-,l+,2+^ 



(4.63) 
(4.64) 

(4.65) 
(4.66) 



The fact that the symmetries ( [4.59| )-( [4.66| ) have to be respected provides a non-trivial 
check on the amplitudes. Since our computation is done by obtaining separate, in general 
non-symmetric, pieces of the amplitude {e.g., the coefficient of a single log in a particular 
channel), it is only after putting them all together that the symmetry becomes manifest. 
Therefore, it is the combination of many non-symmetric ingredients that gives rise to a 
symmetric final answer. An error that spoils the symmetry can be detected by this check. 
We have checked that our amplitudes obey all the required symmetries. 



4.6.3 Numerical comparison 



Ref. |21|] computed the virtual cross section for the H qqQQ process to next-to-leading 
order accuracy (one-loop diagrams, not counting the top quark loop vertex) using a semi- 
numerical approach. They also obtained analytic expressions for the aforementioned cross 
section summed over colors and helicities. 

Using our color- and helicity-decomposed 0-amplitudes presented in our paper, we have 
constructed the same quantity and have compared with their analytical results numerically. 
The result at the phase-space point (|43^ ) was given in table ^, but we have also found 
agreement with their analytical formulae for all the randomly-generated phase-space points 
that we examined. 



5. Conclusions 



In this paper we have presented analytic results for the one-loop amplitudes for a Higgs plus 
four quarks, and for a Higgs plus two quarks and two opposite-helicity gluons. We have 
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obtained the cut-containing and rational pieces of the answer separately, using unitarity 
for the former and on-shell recursion for the latter. We have also shown in specific exam- 
ples how to compute the various ingredients, and presented a way to deal with spurious 
poles without fully eliminating them from the completed-cut terms. Our expressions are 
relatively compact and in agreement with various consistency checks as well as previous 
results. We believe that they will provide a useful input for faster numerical programs 
computing NLO cross sections relevant for the LHC, and will be an important ingredient 
for future higher-point calculations. Together with the NMHV Hqqgg case, and the re- 



maining helicity amplitudes for Hgggg (beyond those already been computed |5l|, |52| 



they provide the one-loop corrections to Higgs plus four partons, and can be used to com- 
pute the gluon fusion contribution to the H + 2 jets final state at the LHC, as well as its 
interference with the vector boson fusion channel. Further NLO studies will be important 
for understanding SM backgrounds, and enhancing the potential for the discovery of new 
physics in the upcoming experiments at the LHC. 
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A. Tree amplitudes 



In this appendix, we record various tree amplitudes entering the main computations and 
results. 



As mentioned in section |3.4.1| , the (j)qq and (jj^qq amplitudes for massless quarks vanish 
by fermion chirality conservation and angular momentum conservation. The (pgg and (t)^ gg 
tree amplitudes are given by 



4^ 



(0), 



A 



(0), 



4°) (<A,1±,2T) = 0, 



A 



-i[12]^ 



(A.l) 
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The 4>qqg and (p^qqg amplitudes are 

4°)(^,ir,2+ 3+) = 0, 



A 



(0), 



1 - o+ o- 



— Z 



(13)^ 
(12) ' 



4°\0t^l-^2+,3-) = O, 
4°)(0t^l_,2+,3+) = 
while the (f)ggg amplitudes are 



[2 3]^ 

w 



4°)(^,i+ 


2+ 


3+) 


4°)(0,i- 


2" 


3+) 


4^0,1- 


2~ 


3-) 


4^^^l- 


2" 


3-) 


4°)(0t,i+ 


2+ 


3-) 


4°)(0t,i+ 


2+ 


3+) 



3+) = 4°)(</.,l-,2+,3+)=0, 

. (12)' 



(2 3) (31) ' 
'[12] [23] [31]' 

if (0,1+, 2-, 3-; 

J12]^ 
[2 3] [31]' 

(12) (2 3) (31) • 



0, 



(A.2) 



(A.3) 



The (pqqQQ and (p^qqQQ tree amplitudes are given by 

|(0) 



4°)f<^,l„2+,3t,4-) = -. ^^^^ 



<^(</'M5,2+ 



4°)(0t,l-,2+ 
For the case of (f)qqgg and 4'^qqgg 



(0), 



A 



3g 

we 



2+ 
2+ 



2+ 
2+ 

9 



2+ 



have 

3+,4+) = 0, 



4;;) 



(12) (34)' 
4°)(</.,ir,2+ 4+,33) 
. [2 3]^ 



[12] [3 4]' 
4+) = -4°)(0t,i-,2+ 4+,3q). 



3 

3- 

3+ 

3- 



■ (14)^2 4) 
(12) (23) (34)' 

(12) (34) (41)' 



4-) = 0, 

4-) = - 

4+) = i 



. [2 3]^ [13] 
112] [34] [41]' 
[2Af 



(A.4) 



(A.5) 



[12] [2 3] [3 4] ' 
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and for the NMHV cases, 



A 



(0), 



1 r ?+ ?, A~ 



4'')(0t,l-,2+ 3+,4 



. (3|(1+4)|2]^(41) 

[2 4] S412 



. (4|(1 + 3)|2]^(13) , . 



1 1 

Sl2 S41 



+ i 



[2 3] S12S123 
^ (1|(2 + 3)|4]^[2 3] 

(13)S123 

(1|(2 + 4)|3]2[24] 



(1|(3 + 4)|2]^ 
(12) [2 4] [2 3] [34] 



1 1 

S12 S23 



(1|(3 + 4)|2]^ 



(14)si2S4i2 [12] (13) (14) (34) ■ 



(A.6) 



(A.7) 



The ampHtudes for (pqgqg and (p^qgqg are, 



A 



(0), 



1 ~ 9+ 3+ 4" 



'1 -^g 



(0), 



1 r 2^ 3+ 4^ 



'1 -^g 



4°)(</.t,lr,2+,3+,4- 



4°)('/'^lg-,2-,3+,4+) 



(14)^ 



(12) (2 3) 



(12)^ 



(34) (41) 
[2 3] 2 



[34] [41] 

. [34]^ 
'[12] [2 3] 



and for the NMHV cases. 



(Ai 



A 



(0), 



A 



(0), 



^ ir,2-,3+,4- 



l„2+,3+ 4+) 



^ (4|(1 + 2)|3]^ (2|(1 + 4)|3]^ 

[12][2 3]si23 [34] [41] S341 ' 
(1|(2 + 3)|4]2 , . (1|(3 + 4)|2]2 



+ 1- 



^q,~ , q,^ , ^(12) (23) 5123 ' "(34) (41) S341 



(A.9) 
(A.IO) 



Our 4>qqgg and (pqgqg tree amphtudes, after dividing by i, agree with the ones (imphcit) 
in refs. [^, 19, 37 1. For the (p^qqgg and 4>^qgqg tree amplitudes, after dividing by i, our 



amplitudes agree with the ones in refs. [18, |l^, but have the opposite sign from those in 
ref. |6^] . (The relative sign between (/) and (p^ amphtudes matters in reconstructing the H 
amplitudes.) Our (pqqg . . . g and (ji^qqg ■ ■ ■ g tree amphtudes are uniformly opposite in sign 



to ref. [|4| 



B. Lj and Ls 1 functions and scalar integrals 



The definitions for the Lq, Li, L2 and Ls_i, Ls^™*^ functions used in the formulae presented 
above are 

Inr lnr + 1-r Inr - (r - l/r)/2 

Lol?-) = , Li(r) = — — , L2(r) = , (B.l) 

1 — r (1 — ry (1 — 



- 51 - 



Ls_i (s, t; rr?) = hh {l - + Lis f 1 - ^ 1 + In 



Ls2r(.,f;m?,mi)= - Lis (l - v) 



- Lio 1 



In 



(B.2) 



m. 



2 1- 



- Lio 1 - 



Li2 ( 1 - Y ) + Li2 ( 1 



2 2 

st 



w - 

2 V 



with 



Li2(x) 



dz 



Ind - 



/o z 

The one-mass and easy two-mass box integrals are related to the Ls_i functions by, 



(B.3) 
(B.4) 



Xr(s,t;m^) 



-2icr 
st 



1 



2 \ e 



- Ls_i(s,t;m ) > , 



-2icr 



st — m\m\ \ e 



1 



2 \ e 



2 \ 



At 



2 \ e 



—mi 



(B.5) 

2 \ en 



— ms 



Lslr(5,i;^!,rn^) 



The one-mass and two-mass triangle integrals are related to each other, 



(B.6) 



— icr 1 / /U 



,2 \ e 



(— s) \ — s 
— icr 1 



2 \ e 



2 \ e 



-S2 



(B.7) 
(B.8) 



and contain terms of the form ln(— Sj)/e; hence their coefficients are dictated by the known 
infrared poles of the amplitude. 



The bubble integral is given by, 

Ms) 



2 \ e 



e(l - 2e) V-s 
and contains a single logarithm, ln(— s), at order e°. 



(B.9) 
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